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Area  of  frame  cross  section 


1  +  4 


pR  2 
+  2  —  X 
Eh  m 


-1 


Frame  width 
Frame  depth 
Young's  modulus 

Distance  between  median  surface  of  shell  and 
cervfcroid  of  frame 

Axial,  circumferential,  and  shear  strains 
Stress  function 
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A  small -deflection  aria  lysis  is  developed  for  the  elastic  i’'.r:erbay 
buckling  of  ring  -  stiffened  cylindrical  shells  in  which  the  influence  of  the 
rings  on  deformations  before  and  during  buckling  is  considered. 

Tests  were  carried  out  with  a  machined,  ring -stiffened  cylinder 
(ER-4B)  subjected  to  external  hydrostatic  pressure.  Collapse  was  initi¬ 
ated  by  elastic  asymmetric  buckling  of  the  shell.  Strain  measurements 
taken  during  the  test  demonstrated  that  the  Southwell  method  of  determin¬ 
ing  buckling  strength  nondestructively  is  applicable  in  the  case  of  inter - 
bay  buckling. 

The  results  of  this  and  three  other  tests  of  machined  cylinders  are 
in  good  .agreement  with  the  theory.  While  the  Von  Mises  theory  is  inac¬ 
curate  for  closely  spaced  rings,  its  continued  application  for  estimating 
elastic  buckling  strength  is  probably  justified  since  it  is  always  conser¬ 
vative  and  can  be  represented  in  a  simple  form. 


INTRODUCTION 

✓ 

One  problem  in  the  field  of  pressure-vessel  design  that  has  been  of 
particular  interest  for  many  years  is  the  elastic  instability  of  thin  cylin¬ 
drical  shells  under  hydrostatic  pressure.  In  practice,  it  is  the  buckling 
which  occurs  between  closely  spaced  ring  stiffeners  as  shown  in  Figure  1. 
Despite  the  considerable  study  which  this-probtem  has  received,  some 
additional  investigation  appeared  needed  in  at  least  two  areas.  One  of 
these  concerned  experimental  work  with  short  shells  since  a  thorough 
evaluation  of  theory  was  still  lacking.  The  other  involved  the  analytical 
study  of  the  influence  of  the  rings  oaThe  buckling  strength.  Several  in- 
vestigatorn  have  considered  this  effect,  but  in  each  case  it  seemed  that 
the  treatment  has  been  either  approximate  or  incomplete. 

As  part  of  a  continuing  program  at  the  Model  Basin  to  study  the 
phenomenon  of  shell  instability,  effort  has  been  directed  toward  eliminat¬ 
ing  some  of  these  deficiencies.  The  present  report,  which  gives  the  re¬ 
sults  of  this  effort,  is  divided  into  three  sections. 


Figure  1  -  Shell  Instability 


The  first  deals  with  the  devel¬ 
opment  of  a  small  deflection  analysis 
which  accounts  for  the  restraint  pro¬ 
vided  by  the  rings.  This  is  done  by 
expressing  the  deflection  both  before 
and  during  buckling  in  trigonometric 
series  form  and  obtaining  a  solution 
through  the  use  of  energy  methods. 
The  special  cases  of  simple  and 
clamped  support  are  also  treated. 

In  addition,  the  significance  of  in¬ 
cluding  the  energy  associated  with 
bending  stresses  developed  prior  to 
buckling  is  examined.  Much  of  the 
work  is  basac|  on  relations  derived 

The  second  section  describes 
tests  of  a  machined  cylinder  desig¬ 
nated  Model  BR-4B.  This  includes 
a  discussion  of  strain  data,  buckling 
strength,  and  mode  of  collapse.  The 
applicability  of  the  Southwell  method 
for  determining  shell  buckling  strength 
is  also  examined. 


In  the  third  section,  the  present  analysis  and  others  are  evaluated 
on  the  basis  of  the  test  results  of  Model  BR-4B  and  of  recent  tests  of 
three  other  machined  cylinders. 


BACKGROUND 


In  view  of  all  the  effort  that  has  been  devoted  to  the  study  of  the 
elastic  instability  of  cylindrical  shells,  it  i9  rather  curious  that  there  has 
yet  to  be  a  thorough  confirmation  of  theory  by  experiment,  at  least  where 
closely  spaced  stiffeners  are  concerned.  Surprisingly  enough,  this  has 


not  been  the  case  with  the  problem 
of  general  instability  wherein  both 
rings  and  shell  undergo  extensive 
deformations  (Figure  Z)  .  It  is  due 
largely  to  the  theoretical  work  of 
Kendrick  1  ’  s:'  3  among  others,  and 
a  rather  extensive  experimental 
program  conducted  at  the  Model 
Basin4'5*6  that  the  elastic  general 
instability  problem  has  for  all  prac¬ 
tical  purposes  been  solved.  This 
progress  is  even  more  striking  when 
it  is  realized  that  no  rigorous  ana¬ 
lytical  solution  was  undertaken 
prior  to  the  work  of  Salerno  and 
Levine7  in  1951, 

'  Wltfct  then,  accounts  for  the 
slower  progress  on  a  problem  which 
would  appear  the  simpler  of  the  two 
and  which  has  been  studied  over  a 
much  longer  period?  Basically  it 
is  the  disparity  which  has  persisted 
in  varying  degree  between  the  phys¬ 
ical  conditions  prevailing  in  the  ex¬ 
periments  and  those  which  have  received  theoretical  consideration. 

First,  let  it  be  said  that  theoretical  development  has  lacked  neither 
variety  nor  distinction.  The  case  of  a  simply  supported  shell  of  finite 
length  under  radial  pressure  was  first  treated  by  Southwell  in  1913.  8 
The  next  year  Von  Mises  presented  a  more  exact  analysis,9  which  he  ex¬ 
tended  in  1929  to  include  end  pressure.  10  In  that  same  year,  Tokugawa11 
published  almost  identical  results  to  those  of  Von  Mises.  While  rightly 
regarded  as  classics  today,  these  analyses  are  not  completely  rigorous 


Figure  2  -  General  Instability 


J  References  are  listed  on  page  101. 
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ill  that  they  do  not  account  for  the  influence  of  the  boundary  conditions  on 
deflections  prior  to  buckling.  This  influence  is  unimportant  for  relatively 
long  shells,  but  it  may  not  be  for  short  shells  such  as  are  embodied  in 
submarine  hulls.  Recognizing  this,  Von  Sanden  and  Tolke12  in  1932  pub¬ 
lished  a  oonipreiieiisi ve  paper  un  ihe  buckling  ox  sheiis  wherein  tbe  eliect 
was  considered.  Their  results  included  not  onJy  the:  case  of  simple  sup¬ 
port  but  that  for  finite  stiffening  rings  as  well.  In  this  latter  case,  how¬ 
ever,  they  assumed  that  during  buckling,  the  rings  merely  provide  simple 
support.  Using  a  different  mathematical  approach,  Kendrick2  also  ob¬ 
tained  a  solution  to  this  problem.  While  neither  analysis  included  the 
rotational  restraint  which  the  rings  presumably  would  provide  during 
buckling,  the  results  showed  that  buckling  pressures  for  short  ring  spac- 
ings  can  be  significantly  higher  than  those  given  by  either  Von  Mises  or 
Tokugawa. 

Salerno  and  Levine ,3' 14  were  apparently  the  first  to  include  the  in¬ 
fluence. of  stiffeners  on  the  buckling  deformations.  Unfortunately,  in 
computing  the  buckling  pressure,  they  neglected  the  prebuckling  defor¬ 
mations  altogether.  Because  of  this  shortcoming  and  because  of  certain 
errors  in  their  energy  expressions,  the  analysis  is  not  considered  cor¬ 
rect.  Others,  notably  Sturm19  and  Nash,16  investigated  the  consequences 
of  having  stiffeners  which  provide  full  fixity.  The  resulting  buckling  pres¬ 
sures  were  much  higher  for  closely  spaced  stiffeners  than  for  the  case  of 
simple  support,  even  though  the  effect  of  the  boundary  conditions  on  the 
prebuckling  deformations  was  neglected. 

Still  others  have  studied  the  possibility  of  "snap-through"  buckling 
at  pressures  much  lower  than  the  small  deflection  analyses  juet  mentioned 
might  predict.  Notable  contributions  to  this  development  of  lax-ge  deflec¬ 
tion  theory  include  the  work  of  Connell, 17' 18  Langhaar  and  Boresi,19  and 
Kempner  and  Crouzet-Pascal.20  These  studies  have  shown  generally  that 
the  phenomenon  of  "snap-through"  in  cylindrical  shells  under  hydrostatic 
pressure  is  possible  and  that  the  associated  buckling  pressure  is  some¬ 
what  influenced  by  geometrical  imperfections.  Experimental  studies  by 
Kir  stein  and  Wenk21  have  borne  out  these  conclusions.  So  far,  however, 
it  has  been  difficult  to  assess  quantitatively  the  imperfections  present  in 
test  structures.  It  has  also  been  recognized  that  such  imperfections  can 
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have  an  important  influence  on  the  stresses  existing  in  the  shell  prior  to 
buckling,  Sturm,15  Bodner  and  Berks,22  and  Galletly  and  Bart,23  for  ex¬ 
ample,  have  shown  that  as  a  result  of  irregularities  in  circularity,  the 
stresses  can  be  so  greatly  increased  that  inelastic  collapse  can  occur  at 
a  pressure  well  below  that  which  elastic  buckling  theory  wouid  predict. 

Inelastic  buckling,  which  is  frequently  the  mode  of  collapse  for  effi¬ 
ciently  designed  shells,  has  also  received  some  attention  in  recent  years. 
For  example,  the  work  of  Gerard, 24,25  Lunchick,26  Nott,27  and  Reynolds  28 
have  shown  that  inelastic  buckling  strength  can  often  be  predicted  with 
reasonable  accuracy  even  though  the  complexity  of  the  problem  virtually 
precludes  a  rigorous  analy?  These  investigators,  it  should  be  noted, 
have  had  the  benefit  of  much  more  experimental  data  than  has  been  avail¬ 
able  for  the  study  of  elastic  buckling. 

It  is  because  of  this  lack  of  data  that  a  thorough  evaluation  of  elas¬ 
tic  buckling  theory  has.  not  been  possible.  Most  of  the  elastic  data  from 
_  Windenburg'a  studies,  29  for  example,  have  beeh  in  the  long- shell  range, 
whereas  data  for  short  shells,  i.e.,  those  with  closely  Bpaced  stiffeners, 
are  extremely  limited,  This  results  from  the  fact  that  the  source  of  much 
of  the  available  data  for  short  shells  has  .been  proof  tests  of  pressure  ves¬ 
sels  designed  for  structural  efficiency  rather  than  for  the  study  of  elastic 
buckling.  Such  structures  are  so  designed  that  collapse  rarely  occurs 
before  stresses  have  exceeded  the  elastic  limit  of  the  material.  Conse¬ 
quently,  tests  of  this  nature  seldom  provide  pertinent  data  for  the  evalua¬ 
tion  of  elastic  buckling  theory,  although  they  have  sometimes  been  used 
for  that  purpose.  On  the  other  hand,  some  cylinders  specifically  designed 
to  buckle  elastically  have  apparently  suffered  premature  failure  by  reason 
of  geometric  imperfections,  residual  fabrication  stresses,  and  other  un¬ 
controllable  factors.  Tests  of  Models  BR-l30  and  BR-531  conducted  at  the 
Model  Basin  are  examples  in  which  this  problem  was  encountered.  How¬ 
ever,  even  if  all  these  experimental  difficulties  had  not  existed,  the  vari¬ 
ous  theories  then  available  still  have  not  considered  realistically  the 
actual  conditions  existing  at  the  shell  boundaries. 

Nearly  all  the  available  data  have  come  from  tests  of  ring-stiffened 
cylinders,  yet  most  analytical  solutions  are  based  on  arbitrary  assump¬ 
tions  regarding  the  boundary  conditions  which  the  stiffeners  impose  on 
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the  shell.  One  remedy  would  be  to  attempt  to  duplicate  in  the  laboratory 
a  set  of  idealized  conditions  as  assumed  for  the  theory.  However,  the 
best  one  could  do  with  such  a  procedure  would  be  to  establish  the  applica¬ 
bility  of  a  certain  formula  in  the  special  case,  for  instance,  of  simple 
support.  The  question  of  what  formula  is  reliable  for  the  case  actually 
encountered  in  practice  would  still  be  unanswered. 

The  recent  successes  in  the  studies  of  general  instability  using 
machined  cylinders  of  high-strength  steel4'5,b  strongly  indicate  that  the 
problems  associated  with  imperfections,  insufficient  yield  strength,  etc., 
could  be  overcome  and  that  elastic  instability  in  shorty  ring -stiffened 
shells  would  be  experimentally  possible.  It  seemed  reasonable,  there¬ 
fore,  to  undertake  an  experimental  study  using  a  machined,  ring-stiffened 
cylinder  designed  to  collapse  by  elastic  instability.  At  the  same  time,  it 
appeared  worthwhile  to  explore  the  possibilities  of  an  analytical  approach 
whereby  the  boundary  conditions  imposed  by  the  rings  could  be  treated 
more  realistically. 


PART  I  -  THEORY 


GENERAL  ANALYSIS 

The  structure  under  consideration  is  a  circular  cylindrical  shell  of 
infinite  length  reinforced  by  regularly  spaced  uniform  ring  frames,  as 
shown  in  Figure  3.  The  shell  is  assumed  to  be  isotropic,  and  its  thick¬ 
ness  is  small  compared  to  its  radius  so  that  the  problem  is  restricted  to 
two  dimensions. 

There  is  ample  experimental  evidence  (e.g.,  Figure  1)  that  under 
hydrostatic  pressure,  asymmetric  (lobar)  buckling  occurs  in  the  mode 
shown  schematically  in  Figure  4.  Such  a  pattern,  it  will  be  seen,  repeats 
itself  at  every  other  frame.  Consequently,  in  computing  the  total  poten¬ 
tial  of  the  cylinder  in  its  buckled  state,  it  is  convenient  to  carry  out  the 
calculation  over  any  two  adjacent  frame  spaces. 

The  total  potential  U,p  of  the  system  is  defined  by 

UT  =  us  +  u£  -  W  f  1  ] 


6 


Figure  3  -  Stiffened  Cylinder  and 
Coordinate  System 


Figure  4  -  Asymmetric  Shell 
Buckling  Mode 


where  Ug  is  the  strain  energy  of  the  shell, 

Uj  is  the  strain  energy  of  the  frames,  and 

W  is  the  work  done  by  the  external  pressure. 

Referring  to  Appendixes  A,  B,  and  C,  where  these  three  quantities  are 
derived,  we  have  for  the  total  potential 
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and  all  terms  beyond  the  third  order  in  the  displacements  have  been  dis¬ 
carded,  The  quantities  appearing  in  Equation  [2]  are  dofined  as  follows: 

E  is  Young's  Modulus, 

h  is  the  shell  thickness, 

R  is  the  radius  to  the  shell  middle  surface, 

Lf  is  the  frame  spacing, 

v  is  Poisson's  ratio, 

Aj  is  the  area  of  the  frame  cross  section, 

*xG  is  ^ie  moment  of  inertia  of  a  frame  about  the  centroid  in 
its  plane  of  curvature  (Figure  3), 
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Izq  is  the  moment  of  inertia  of  a  frame  about  the  centroid  out 
of  its  plane  of  curvature  (Figure  3), 

e  is  the  distance  from  the  frame  centroid  to  the  neutral  axis 
of  shell  (Figure  3k 


K  is  the  torsion  constant  of  the  frame  (Appendix  B), 
p  is  the  hydrostatic  pressure  (positive  outward), 


0  and  x  are  the  angular  and  axial  coordinates  (x  positive  to  the 

iv-  '  5  i  i-  *  •  !  5  ik.  /.H'tl 

right),  and 

u,  v  and  w  are  the  axial,  tangential  and  radial  displacements,  respec¬ 
tively  (u  positive  to  the  right,  w  positive  outward). 


The  subscripts  indicate  differentiation. 


During  buckling,  the  system  passes  from  an  initial  or  prebuckling 
equilibrium  state,  in  which  all  deformations  are  axi symmetric,  to  the 
buckled  Or  asymmetric  state.  The  change  in  total  potential  accompanying 


this  process  may  be  called  AUrp  .  The  final  buckled  state  will  be  described 
by  the  deflections  ufl  Vp  and  Wp  which  are  the  total  of  the  displacements 
developed  from  the  initial  application  of  pressure  and  are  given  by 


uF  =  u(x)  +  u  ( x,  0 ) 

Vp  =  v(x,0)  [3] 

Wp  =  w(x)  +  w  (x,  0) 


where  u,  v,  w  are  the  buckling  displacements  and  u,  w  are  the  initial  or 
prebuckling  displacements,  being  axisymmetric  in  form  for  the  case  of 
hydrostatic  pressure,  and  satisfying  the  conditions 
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To  find  the  change  in  total  potential,  we  replace  u,  v,  w  in  Equation  [2] 
by  Up,  Vp,  Wp  and  subtract  the  total  potential  for  the  initial  equilibrium 
state.  With  only  linear  and  quadratic  terms  retained,  the  result  is 
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where  M  and  N  now  involve  only  the  buckling  displacementfl  v  and  w. 

A 8  explained  in  Appendix  D,  u  and  w  can  be  expressed  by  the  closed 
form  solution  of  Von  Sanden  and  Gunther32  for  the  linear  case,  or  by  the 
more  exact  nonlinear  solution  of  Pulos  and  Salerno,33  '  However,  as  can 
be  seen  from  Appendix  D,  manipulation  of  either  of  these  solutions  could 
be  extremely  cumbersome.  Kendrick2  provides  convincing  evidence  of 
this  in  an  analysis  of  the  same  buckling  problem, 

A  simpler  approach  is  to  express  u  and  w  in  terms  of  infinite  trig¬ 
onometric  series,  as  indicated  by  Von  Sanden  and  Tolke.12  This  has 
been  carried  out  in  Appendix  ,D  with  the  following  results: 
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continued 


where,  for  the  equivalent  of  the  Salerno  and  Pulos  solution, 
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If  the  pressure  term  in  the  denominator  of  am  is  neglected,  we  have  the 
equivalent  of  the  Von  Sanden  and  Gunther  solution:* 
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This  approximate  form  is  found  to  be  sufficiently  accurate  for  the  range 
of  geometry  to  be  considered,  particularly  since  this  analysis  is  not  con¬ 
cerned  with  an  examination  of  stresses  and  deflections  at  discrete  points. 

00 

Appendix  D  also  gives  an  alternative  means  of  computing  £  am  for 
the  approximate  form:  rn~-a> 


J3  /  sinh  P  +  sin  P  \ 
2  \  cosh  (3  -  cosp  / 


19] 


small  discrepancy  between  the  series  and  closed  form  solutions  re¬ 
sults  from  neglecting  the  faying  width  of  the  frame,  as  explained  in 
Appendix  D. 


13 


Except  for  very  large  values  of  p,  however,  the  convergence  is  so  rapid 
that  the  series  form  is  often  the  more  convenient. 

Having  the  necessary  Equations  [5]  and  [61.  it  is  possible  at  this 
point  to  apply  the  principle  of  stationary  potential  and,  through  the  meth¬ 
ods  of  variational  calculus,  proceed  to  a  solution  of  the  problem.  The 
most  rigorous  procedure  would  be  to  formulate  the  partial  differential 
equations  of  equilibrium  and  then  attempt  to  solve  them.  This  approach 
is  rejected  here  primarily  because  of  the  lengthy  task  which  would  be  in¬ 
volved  in  obtaining  the  differential  equations.  The  work  of  Von  Sanden 
and  Tdlke12  indicates  that  an  exact  solution  to  the  differential  equations 
may  be  entirely  possible,  provided  the  initial  deflections  are  expressed 
in  their  trigonometric  series  form.  A  simpler  procedure  which  will  be 
followed  here  is  to  apply  the  Ritz  method  whereby  displacement  functions 
having  arbitrary  coefficients  are  assumed  and  the  problem  of  satisfying 
the  condition  for  stationary  potential  ie  reduced  to  the  solution  of  a  system 
of  algebraic  equations.  ,  ,v 

Since  the  success  of  this  method  depends  upon  the  degree  to  which 
the  assumed  functions  approximate  the  exact  buckling  shape,  it  is  impor¬ 
tant  to  choose  a  set  of  functions  which  permit  a  wide  variation  in  shape. 
Accordingly.  th«  dinnlar  Aments  are  annnmnd  an  fnllnwn; 

to 

u  =  cosne  £  um  008  (“r“) 

m=  1 


v  =  sin  n0  jr  vm 
m=] 


where  m  and  n  are  positive  integers  and  u^  vm,  wm  are  the  arbitrary 
coefficients. 

It  can  be  seen  from  Figure  5  that  with  appropriate  adjustment  of 
the  coefficients,  buckling  configurations  satisfying  simple  support,  r 
tial  fixity,  or  full  fixity  can  be  generated.  It  should  be  recognize 
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B 


Figure  5  -  Buckling  Configurations 

frrr  TKraa  Haora* n  rtf  Pi  yHv 


although  v  and  w  vanish  at  each 
frame  for  any  set  of  coefficients,  no 
generality  is  lost,  since  this  is  one 
of  the  conditions  which  define  the 
present  buckling  problem.  Buckling 
configurations  involving  radial  and 
tangential  deflections  of  the  rings 
properly  belong  to  the  category  of 
overall  buckling  or  general  instabil¬ 
ity  which  is  outside  the  scope  of  this 
analysis. 

While  the  procedure  for  solving 
the  problem  in  this  general  form  is 
straightforward,  considerable  effort 
would  be  involved,  and  the  results 
would  probably  be  unmanageable, 
it  is  desirable,  therefore,  to  intro¬ 
duce  simplifying  approximations 
where  possible  before  proceeding 
further.  The  sections  that  follow 
cousiut.!  w^iy  tivrco  c «x t* c s  of 


interest:  partial  fixity  (finite  rings), 
simple  support  and  complete  fixity  (infinite  rings).  Certain  approxima¬ 
tions  and  methods  peculiar  to  each  case  are  developed  which  lead  to  sim¬ 
plified  results. 


RINGS  OF  FINITE  RIGIDITY 

In  this  case,  rings  of  practical  size  are  considered.  The  basic 
assumption  is  that  the  various  ring  properties  are  of  the  same  order  of 
magnitude  as  comparable  quantities  for  the  shell.  The  basic  approxima¬ 
tion  is  to  eliminate  the  unknown  deflections  u  and  v  through  the  use  of  a 
stress  function,  as  explained  in  Appendix  E.  The  complexity  of  the  buck¬ 
ling  equations  is  thereby  reduced  by  a  factor  of  3.  Since  this  procedure 
constrains  «  and  v,  it  leads  inherently  to  a  higher  buckling  pressure  than 
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f 


would  result  if  complete  generality  were  preserved.  However,  other  in¬ 
vestigators  have  found  that  the  error  thus  introduced  is  extremely  small 
for  the  case  of  a  simply  supported  shell  under  hydrostatic  loading.  The 
assumption  in  this  analysis  is  that  the  error  will  also  be  small  in  the  case 
of  ring  support  where  the  rings  are  of  practical  size  .  In  the  case  of  a  fully 
clamped  shell  (corresponding  to  a  ring  of  infinite  rigidity),  the  assump¬ 
tion  is  not  valid  for  the  assumed  buckling  configuration.  Equation  [  10], 
as  will  be  discussed  in  a  later  section. 

From  Appendix  E,  the  equations  defining  the  stress  function  F(x, 0) 
and  the  deflections  are 
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The  remaining  unknown  deflection  w  is  still  assumed  to  have  the  trigono¬ 
metric  series  form  (Equation  [  10]).  If  this  series  is  now  substituted  in 
Equations  [11],  the  result  is 
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When  these  quantities  are  introduced  in  Equation  [5],  certain  simplifica¬ 
tions  arise.  It  can  be  seen  that 
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when  substitutions  are  made  from  Equation  [11].  Furthermore, 
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because  of  the  periodicity  of  F  in  both  the  x  and  9  directions,  as  indicated 
by  Equation  [12].  Hence  Equation  [13]  becomes 
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Equations  [  12]  also  provide  a  useful  approximation  which  is  valid  for 
short  shells.  The  quantity  N,  which  appears  repeatedly  in  Equation  [5], 
can  be  written: 
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Since  <j>m  S  1  and,  for  short  shells,  X^,  »  1,  little  accuracy  will  be  lost 
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if  one  makes  the  approximation1' 
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In  addition,  it  will  be  seen  that  several  terms  in  the  frame  energy  inte¬ 
grals  vanish  at  the  frames  according  to  Equations  [10]  and  [12]. 
Accordingly,  the  simplified  form  of  Equation  [5]  is 
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*The  accuracy  of  this  approximation  depends  not  on  the  size  of  n  but 
only  on  the  size  of  Xm.  Its  range  of  validity,  therefore,  includes  the  case 
of  axi symmetric  buckling  of  a  short  shell  (n=  0,  =  1)  hut  not  that  of 

asymmetric  buckling  of  a  long  shell  where  Xm  may  not  be  large. 
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where  AUjj  is  the  portion  of  the  bending  strain  energy  which  involves  the 
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initial  deflections.  It  is  given  by 
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will  vanish  on  integration. 


19 


AUb 


2w: 


X6  R2 


ee 


+  2vwxxw00 


-  — - . 2  J  f  f  {t[<3-5v> 

24R(l-v2)  J0  J0  & 

+  2wx  [r2uxxwxx  +  wu^  WQe  -(wx-^)(^+2*Rwxx) 


] 


+  (1  -  v)wx0  (u*0  4  4  Rw^  2R  (2^ 


„  .  we  vxx  . 

wxx  +  — - +  uxxw 


.) 


T  PV  ,B)}]  +  "*  [zr2*5x  +  U'SWwJe  +  &*MWB] 


+  2R5„rRw_w....,^?l\ 

—  L  '*  K  J  I 


[19] 


dx  d0 


After  the  eerie*  expressions,  Equation*  [6]  and  [  12],  are  intro¬ 
duced,  Equation  [18]  can  be  integrated.  In  *o  doing,  it  should  be  noted 
that  because  the  function*  of  Equation  [  12]  are  orthogonal  in  the  interval 
0  5  x  S  2Lf ,  coupling  of  different  buckling  mode*  (i ,  j)  will  arise  only 
from  the  frame  energy  and  from  shell  energy  terms  involving  product* 
of  the  buckling  modes  with  the  initial  deflection*.  The  resulting  inte¬ 
grated  form  of  Equation  [  18]  is 
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i-1  j-1 


“  “Ki+j)  ‘ft2  +  2Xj  (Xi  +  Xj  )fij  ’  qijJ  }  +  AUt 


where: 

«y  =  l+(-l)‘+J 

rm  1  1  +  4,m  +  [*m(> -'t>mXl  + v)‘  -  (1+ 2v)] 

Am 

^ij  =  d  -  fy)  +  <j>j  ( 1  “  4>j) 

qtj  =  4-  4».(1  -<{>.)[!  -  {1  -  v)(4».  -  «h)  -  ♦j+jd-  3v)] 


AUfc  has  not  been  given  explicitly  at  this  point  Bince  it  will  be  necessary  to 
examine  only  a  few  of  its  characteristic  terms.  For  example,  we  find 
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dxd0 
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irpLf(2-v)r  h2n2 
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in-  i 


This  can  be  compared  with  other  terms,  designated  by  Z,  appearing  in 

AU<p : 


Z  = 


co  a/ 


m=  1 


2n  wjm  -  (2  -  v)  T), 


i=l  j=l 


n2wiwjl»i(i.j)-»i(i+j)]][22] 


h2n2 


It  is  immediately  seen  that  if  - 5-  «  1  ,  which  is  a  reasonable  approxi- 

6R2 

mation  for  thin  shells,  Ij,  is  negligible  compared  with  Z. 


ZLf  f 2 7T 
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Eh3  r  cut  r 

24R{l-v2)  in  J, 


w  w wn 

2(1  -  v)  —■  ~~  dxd0 
R 


~ 


CO  CD, 


<z  ■ v)  ”  IsH  J  £  £ n2wi  Wj  x* |a|<i+j)  <Xi + Xj 


i=i  j=i 


+  a*(H>(Xi“\)>l 


[23] 


Comparing  l£  with  Z,  we  find  that  Ig  is  negligible  if  — —  (Xj  +  XcJX^ 

h2  6RZ  J 

and  - — -  ( X •  ~  X  * )  X  ■  are  negligible  compared  with  unity.  This  will  hold 

6R2  J 

for  short  thin  shells  except  when  i  or  j  are  large.  In  such  cases,  however, 
it  will  be  seen  from  Equation  [20]  that  ai, . ,  ..  or  ai, ,  will  be  extremely 
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small  so  that  all  terms  in  which  either  appears  as  a  multiple  will  be 
negligible . 

A  third  term  is 


h  = 


Eh‘ 


/•2Lf  /*2tt 

I  vR 

'0  J0 


24R ( 1  - v  ) 


l^v2)  J 


Wjjx  w^dxd0 


i/irpLi 


CO  CO 


2-'M12Ri?b(l-v2)]  E  E™ iWJXiXJ  (ai(i-j)IXi*Xj)2 

[24] 


+  ai(i+.t)(Xi  +  Xj»] 


A  comparison  with  Z  reveals  that  I3  will  also  be  negligible  eo  long  as 

h2X.  X;  0  h2X,X.  9 

- —  (X.  -  X;)”  and  - - — - —  (X.  +  X-)“  are  negligible  compared 

12RZ(  1-v2)  V  .  t  . „  12R2  ( 1-v2)  1  J  ...  .  . 

with  n2.  As  In  the  case  of  I2  ,  when  this  does  not  hold,  the  question  no 

longer  has  importance. 

It  can  be  shown  further  that  each  term  appearing  in  AU^  will  be  one 
of  these  three  general  tvoes.  Hence,  the  entire  integral  AUu  is  negligi- 
Die.  m  otner  woras,  tne  pending  stresses  associated  wan  tne  preouciuing 
deformations  can  be  neglected. 

Inspection  of  Equation  [20]  reveals  some  additional  approximations 
which  can  be  made.  It  will  be  seen  from  Equation  [  12]  that  Um  is  negli¬ 
gible  compared  with  Xm  since  4>m  is  always  less  than  unity  and  Xm  is 
much  greater  than  unity.  Hence, 

wm<Um  +  Xm>*  wm  Xm  lZ5J 


It  also  follows  that 


^mm 


w. 


m 


[26] 
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2  lit 

If  n  .»  1,  we  can  also  neglect  ,  so  that 


ai(i-j)  "  ^Aj,Ai  "  Aj,fij  f  qij5  “  aj{i+j)  ln“  +  2Xj(Ai  +  Xj)  *ij  “  ^ij1 

“  ai( i . j>  t n?'  -  2X j a}( r j )  i >z + 2Xj <vw 


[27] 


When  i  »  j  ,  this  approximation  will  not  hold.  However,  as  explained 
for  1^  and  I^,  the  entire  expression  (Equation  [27])  in  such  cases  is 
negligible. 

We  have  then  for  the  approximate  form  of  Equation  [20]: 
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[28] 


inis  is  generally  a  valid  assumption  for  asymmetric  buckling  of  short 
shells;  for  axisymmetric  buckling,  is  identically  zero. 
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npL{ 


UJ  QJ 

ZZ^i  {ai(i-j)[n2  '  2Xj<xi  -  Wji 
i=1  j"1  [28] 


4<i+j)  ["Z+  ZXjOli+V'ij]} 


continued 


The  buckling  condition  is  now  obtained  using  the  principle  of  station¬ 
ary  potential.  Minimizing  AU^  with  respect  to  the  Wj  and  multiplying  the 

R 

resulting  equations  by  — — — —  ,  one  obtains  the  simultaneous  equations: 

ft  E I .if  h 


8(AUt)  r 
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12R^l-v^)4> 
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i=  1 


ai.;  [n2  +  g..  (\.  -f  X-)2l\ 
2<l-j)  1  J1  1  J 
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,  .  .  pR 

where  dj  is  -  . 

Eh 

u  ;  „  _  „  «  „«««  --  ~  -  -  x*  tt' *. : r  t  n  1  i\.  .n  -  1 1  —  -  .1  -  - 

It  IS  2.ZI  intC-  Cuiai w^u«jul>c  ux  uv|uaLiuiio  [  *-  7  J  axi  muuc » 

corresponding  to  odd  values  of  i  and  j  are  coupled  and  the  same  is  true  of 
all  even  modes,  but.  that  no  coupling  exists  between  odd  and  even  modes. 
This  arises  from  the  two  consistent  definitions 


sji  =  1  + 

=  aa(i+j) 


i+j  _ 


unless  i,j  odd  or  i,j  even  [30] 


The  problem  thus  can  be  grouped  into  two  sets  of  equations  corresponding 
to  the  odd  and  the  even  modes,  respectively.  It  is  clear  that  only  the  odd 
modes  need  be  considered  here  since  the  even  modes  lead  to  buckling  con¬ 
figurations  of  the  type  (shown  in  Figure  6)  which  repeats  itself  at  each 
frame  and-.is  x)f>-no  practical  significance.  When  the  problem  te  reatrleted 
to  the  odd  modes.  Equations  [29]  become: 
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i»  j  =  1,  3,  5 


Figure  6  -  Even  Mode  Buckling 

If  j  is  given  an  upper  limit  J,  Equations  [31]  can  be  written  in  the  form 

■L,  r  J  , 

BJwj+  *Z7biiwi  +  ,MClwi  +  A=iiwi  I  *  0  [3Z] 


£ bji  wi  +  +  [cj  wj  +  wi  ]  * 
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A  requirement  for  the  existence  of  a  nontrivial  solution  to  these 
equations  is  that  the  determinant  formed  from  the  coefficients  of  the  wj 
must  vanish.  This  rrmditinn  ran  he  written  in  matrix  notation  as  follows: 


where: 


{[R]  +  [.Cl  }  [X]  =  0 


[X]  = 


Bl  +  bll’  b13*  bl  5  ’ 


b3i *  B3+b33' 


b51  » 


B5+  b55  ’ 


bJl' 


BJ  +  bJJ 


Ci  +  cir 

C31  *  C 


C51  ' 


Both  [B]  and  [C]  are  symmetrical  matrixes.  The  equation  has  — —  real 
roots,  all  of  which  are  negative,  since  pressure  was  defined  positive  when 
external.  Of  these  roots,  the  one  (4^)  whose  absolute  value  is  least  defines 
Pn,  the  buckling  pressure  for  n  circumferential  lobee.  Since  n  is  an  arbi¬ 
trary  integer,  it  is  necessary  to  obtain  pn  for  several  values  of  n  until  the 
minimum,  or  critical  pressure,  pc  is  determined.  The  degree  of  accuracy 
depends,  of  course,  on  the  magnitude  of  J.  As  J  increases,  pc  should  con¬ 
verge  toward  a  limiting  value. 

While  the  solution  of  Equation  [33]  appears  reasonably  adaptable  to 
a  high-speed  digital  computer,  it  becomes  unmanageable  for  desk  calcula¬ 
tion  if  J  is  much  greater  than  5.  To  overcome  this  difficulty,  an  approxi¬ 
mate  form  has  been  found  whereby  the  computational  task  is  considerably 
reduced.  With  suitable  rearrangement,  Equation  [32]  can  be  put  in  the 
form: 


[Bj  +  i|iCj)Wj  +  Gj 


wiAi  +  Hj 


^wiui  +  *  A.  wi  dji  =  0 

i=i  i=i 


where  Bj  and  Cj  are  as  before  defined  for  Equation  [32]  and 


j  R2(l+|)3Lfh  L1  R 


i-— •  in  \  a.ii  -  -in  -AM-t-n  Ui 


'JIM*'!* 


+ - A  n2 

(l  +  v)R2(l  +  |)3Lfh  1+|  \  j  L\  R2Af  / 

-|(I  +  l)]-n2|  uj} 

=  2  2\Cf  fxj  [l  -  |  (n2  -  1)1  +  n2  U,}  4  't'<2  ~  n2  (u,  -  tOl ) 

R2(l+|)3Lfh  ^  J  L  R  J  U  1  +  |  \J  R/ 


dji  =  ^Vii{a.(i+j)[n2+  IjilAj+Aj)2]-  a.,  [„2+  -  A.)2]} 


To  simplify  the  equations,  we  neglect  all  the  d-  for  which  j  and  i  are 
unequal,  so  that  Equation  [34]  becomes 


Dj  wj  +  gj  JT'Vi + hj  ui 

i~l  i=! 


j.i  =  1.  3,  5 


35] 


wher<s: 


Dj  =  Bj  +  »MCj  +  djj) 


3  (2  -  v)q  r  ,  2  .  x  2\ 

djj  =  - 2 - 


[aj(n2  +  4ljj  X2 )  -  n2] 


It  can  be  seen  from  Equation  [8]  that  the  approximation  is  not  unreason¬ 
able,  Because  decreases  rapidly  with  increasing  m,  most  of  the  am 

will  be  negligible.  In  many  cases,  it  would,  in  fact,  be  neceBsary  to  re¬ 


tain  aQ  only ;  A  more  conclusive  test  of  the  apprbximatidft,  how^vir,  is  W 


comparison  of  numerical  results.  Since  am  approaches  zero  more  slow¬ 
ly  as  p  is  increased,  a  demonstration  that  the  approximation  is  valid  for 
a  large  value  of  p  should  be  sufficient  proof  of  its  validity  for  smaller 


V  ■ 


cal  example. 

9 

With  Ijj  expressed  by  Equation  [20]  and  n  replaced  by 
a  more  concise  expression  for  Dj  is  obtained: 


Dj  =  Bj  +  |  [*j  +  n2{2  -  (2-vH  [1  -  aj(l+  S^2)]}]  [36] 


The  advantage  of  the  form  of  Equations  [35]  over  that  of  Equations  [32] 

can  now  be  demonstrated.  With  the  stipulation  that  \\i  is  such  that  all  D- 

J 

are  different  from  zero,  one  can  multiply  Equation  [34]  by  X./D.  and 

3  j 

carry  out  a  summation  on  j  from  1  to  J: 
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A  second  equation  is  obtained  by  multiplying  Equation  [35]  by  Uj/Dj  and 
performing  the  same  summation: 


A  r  v’Hiuo  v’uiui 
2>iui  L1  + X  —  J+4  — 

i=i  1=1  1=1 


A,  G.  U. 


J 

z 

1=1 


W;X-  =  0 
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38] 


The  unknown  coefficients  wi  are  eliminated  by  combining  Equations  [37] 
and  [38],  and  the  result  is  a  single  equation. 


r  ^  Hiuilr  +  y  Gi Xi  1  y  GiUi  yHiXi  _ 


=  0  [39] 


which  can  be  solved  by  plotting  the  left-hand  side  versus  4*  •  Such  a  plot 

will  V)  vp  A  7.pr  n-int'Prrp  -rvt  c  n-r>  fr\y  nf  Wwa  ortna  ti/vr*  a  r* 

2  ir’~'  ^ - - - 

equal  number  of  asymptotes  corresponding  tc  the  vanishing  of  each  of  the 

D-.  The  root  whose  absolute  value  is  a  minimum  will  lie  between  the  first 

two  asymptotes  (i  =  1,  3).  As  in  the  case  of  Equation  [33],  it  is  necessary 

to  try  several  values  of  n  until  the  critical  pressure  p  is  determined. 

The  accuracy  of  the  result  will,  of  course,  improve  as  J  is  increased. 

In  the  case  where  an  initial  value  of  J  has  been  found  to  be  insufficient, 

one  can  proceed  to  the  next  succeeding  J  (i.e.  ,  increase  the  order  of  the 

equations  by  one)  with  far  less  labor  than  would  be  involved  when  working 

with  the  matrix.  Equation  [  33]  . 


SIMPLE  SUPPORT  CONDITIONS 

The  special  case  of  a  shell  having  regularly  spaced  simple  supports 

.  .  - 

instead  oi  unite  rings  can  be  obtained  quite  easily  from  the  general  equa¬ 
tions.  For  this  case,  the  conditions  are 

*zG  =  *xG  =  K  =  0 

Af  =  CO 

w  =  w  =  0  at  x  =  0,  Lf  [40] 
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[40] 

continued 


’ll 


il  = 


CO  ' 


m 


2  /cosh  P  -  cos  p  | 
p  \  einii  p  +  tun  p; 


if  pressure  term  in  a  can 
be  neglected 


Substitution  of  these  relations  directly  into  Equation  [32]  will  not  lead  to 
the  correct  result  since  some  of  the  teims  in  that  equation  have  been  mul¬ 
tiplied  by  /i£ ,  and  this  would  be  equivalent  to  dividing  by  aero.  Instead 
we  return  to  Equation  [18]  and  note  that  because  of  the  relations  in  Equa¬ 
tion  [40],  ail  frame  integrals  in  which  w  appears  as  a  multiplier  as  well 
as  those  multiplied  by  Izq,  K  and  Ixq  will  now  be  zero.  It  is  then  evi¬ 
dent  that  Equation  [  34]  for  this  case  will  reduce  to 


J 


where  the  Bj,  Cj,  and  djj  are  defined  in  Equations  [32]  and  [34].  This 
equation  can,  of  course,  be  Solved  in  exactly  the  same  manner  as  Equation 


Equation  [3b],  one  arrives  at  a  simple  linear  equation  for  this  case.  That 
is,  if  all  dj£  are  neglected  where  j  and  i  are  unequal,  the  resulting  equa¬ 
tion  has  the  uncoupled  form: 


12RZ(l-vZ)4>? 

_  _  J 


aj  ^  + 


[42] 


This  solution*  gives  a  minimum  buckling  pressure  for  j  =  1 ,  which  is  the 


*It  may  be  of  some  interest  to  compare  Equation  [42]  with  the  equiva¬ 
lent  solution  of  Von  Mi  see  (Equation  [8]  of  Reference  29),  which  in  the 
terminology  of  this  report  is: 
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first  asymptote  of  Equation  [39]  with  set  equal  to  . 

It  is  now  possible  to  find  the  value  of  n  for  which  pn  will  be  a  mini- 

mum .  If  a2  is  replaced  by  X["  j  *s-  l  42]  >  Prt  w*iX  be  effec¬ 

tively  minimised  by  setting  9p  / #4^  equal  to  aero.  The  resulting  equa¬ 
tion  is 


^(Footnote  continued  from  preceding  page) 

r  _  h2xf 


.fr+!+__pi_| 

R  L  12RZ(1  -  vZ)4>Z  J 


^n 


ll  .  2 

—  r  n 


not  affect  the  prebuckling  deformation* .  The  Von  Miles  pressure*  will, 
in  general,  be  lese  than  tfcoee  given  by  Equation  [42  ],  except  where 
aj(l+  tyj )  >  1 .  The  two  solutions  are  identical  fur  p  »  1 ,  where  T)w  ap¬ 
proaches  zero,  but  are  not  accurate  in  this  range.  For  such  a  case,  the 
y««  \/{  ■  •  ■  formula  [  6]  of  Refer  sues  29  should  bo  used. 

A  comparable  solution  ia  that  of  Von  Sanden  and  |Tolkeu  which,  for 
the  case  of  simple  support,  can  be  written 

„  i>2xf  1 


^  = 


RL  12RZ(1  -  v  )  <j>Z  J 


cosh  £  -  cos 
,  4  2  2 
_  =  - - - 

*  P  .  ,  B  3 

_  sinh  ~  +  sin  ■—  __ 

In  deriving  this  formula,  Von  Sanden  and  Tolke  took  into  account  the  ef¬ 
fect  of  the  boundary  conditions  on  the  prebuckling  deformations.  However, 
because  they  did  not  use  the  exact  Fourier  coefficients  for  w  ,  their  equa¬ 
tion  is  not  identical  with  Equation  [42]  of  the  present  analysis. 
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obtained  directly  from  Figure  7,  is  shown  in  Figure  8  together  with  a  com 
parable  curve  obtained  from  Windenburg's  Equation  [24]^  of  Reference 
29,  based  on  the  Von  Mises  solution,  ft  should  be  noted  that  according  to 
Equation  [43],  the  condition  .for  axisymmetric  buckling  { -  I,  or  n  -  0) 
occurs  where  the  value  of  p  is  approximately  1.807.  In  cases  where  p  Is 
less  than  1.807,  Equation  [43]  has  no  physical  significance,  n  being  zero 
for  all  such  cases.  The  plot  of  Equation  [43]  in  Figure  8  shows  that  as  p 
ia  reduced,  n  does  not  increase  indefinitely.  Evidently,  it  attains  some 
limiting  value  after  which  it  decreases  until  the  axisymmetric  condition 
is  reached. 

With  appropriate  values  of  <J>  ^  determined  from  Figure  7,  a  plot  of 
Equation  [44]  was  obtained  and  is  shown  as  the  solid  curve  of  Figure  9 


Figure  9  -  Ruckling  Pressure  as  a  Function  of  p  for 
Simple  Support  (v  =  0.3) 


pr  /R\2 

where  f~)  is  plotted  a?  a  function  of  (3.  Also  shown  is  a  curve  for 
a  xi  symmetric  buckling  representing  the  well-known  equation 


fc/R  \* 

E  \h  / 


/  ? 
Nf  3(1  -  v*"> 


(-H2 + 4 


2p 


[«] 


for  (3  <  7r  ,  This  equation  is  readily  obtained  from  Equation  [42]  it  n  is 
set  equal  to  zero  and  is  the  same  &f>  that  for  the  buckling  of  a  cylinder 
under  end  thrust*4  where  the  length  of  a  full  longitudinal  wave  is  equal  to 
2L,£  .  It  will  be  seen  that  the  two  curves  intersect  at  the  point  (3  s  1.807, 

P  2 

=  -  1.315  .  Since  Equation  [45]  is  not  applicable  for  p  >1.807,  the 
curve  is  shown  dashed  in  that  region.  Similarly,  the  curve  of  Equation 
[44]  xb  shown  dashed  where  P  <1.807.  Figure  9  also  includes  a  plot  of 
the  familiar  equation  of  Windcnburg29  (shown  dotted),  often  referred  to  as 

"EMB  Formula  [  10] . "  This  is  am  approximate  . r.apjre sentatlonqf  theVon 
Mises  solution  (see  earlier  footnote)  in  minimized  form  which,  in  the  ter¬ 
minology  of  this  report,  can  be  expressed  by* 


E  \h  / 


«  -  0.8175 


r*»j 


This  equation  and  all  three  curves  of  Figure  9  are  based  on  a  value  of  v  of 
0.  3  .  Equation  [46]  is  more  conservative  than  Equation  [44]  in  the  low-p 
end  of  the  asymmetric  -buckling  range  because,  as  was  pointed  out  in  a 
previous  footnote,  the  Von  Mises  analysis  neglects  the  effect  of  boundary 
conditions  on  the  initial  (prebuckling)  deformations. 

.It  should  be  mentioned  that  the  "beam-column”  effect  represented  by 
the  factor  p/pm  appearing  in  was  found  to  have  a  negligible  influence 


*The  formula  is  better  known  in  its  original  form: 

5/2 


:'60(|r) 


L  n  ,,/h  \*/2 

JR  "  0,'*5\2R/ 
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on  the  curves  of  Figures  7,  3,  and  9.  The  approximate  form  of  Equation 
[40]  la,  therefore,  sufficiently  accurate,  anti  was  used  in  the  construction 
of  the  curves. 

RINGS  OF  INFINITE  RIGIDITY  (C DAMPED  SUPPORT  CONDITIONS) 

Another  special  case  of  some  interest  is  that,  of  a  shell  having  infi¬ 
nitely  rigid  rings.  As  indicated  previously,  a  solution  for  this  case  cannot 
be  obtained  from  the  equations  for  finite  rings  by  letting  the  frame  rigidi¬ 
ties  become  infinite.  This  is  because  the  use  of  the  stress  function  approx 
imation  has  imposed  certain  constraints  upon  u  and  v  which  can  give  rise 
to  serious  errors  when  the  frame  parameters  become  very  large.  It  is 
preferable  instead  to  proceed  from  the  more  general  Equation  15],  allow¬ 
ing  full  variability  of  u  and  v  (consistent  with  the  boundary  conditions)  as 
well  as  of  w.  Since  the  shell  is  now  fully  clamped  at  regular  intervals 
against  any  deformations,  the  boundary  conditions  are* 


w  =  w  =  v  =  wx  =  Ug  =  0  at  x  ss  0,  hf  1 47  ] 


A  il  .  *  •  % 


tiiss  bu;iu:;iuu  uu  w  i  C4”ifC5 


at-  -  * 
bUAl 


m=-co 


2_  /  cosh  (3  -  cos  ft  \ 

(3  \  sinh  {}  +  sin  p  /  ' 


if  pressure  term  in  a  can  be 
neglected 


The  effect  of  imposing  the  conditions  of  Equation  [47]  is  to  cause  ail 
frame  integrals  in  Equation  [5]  to  vanish.  If  it  is  assumed  as  before  that 
the  bending  energy  AU^  involving  the  initial  deflections  can  be  neglected. 
Equation  [5]  becomes 


'’’These  conditions  permit  axi symmetric  translation  at  the  boundaries. 


:hR  f 2i,J t2" I  2,1 WT vo\  r*v9\  fi-vM  ua\‘  1 

nr,J.  J,  h+(— ) +  2 t)  bd# 


*rr„* 


,  W  +  V 


W  +  V, 


■  \  #  v  V/ 


*2L#  r  ,  ->  +  vg  v2 


Z4R(  1-v  )  J0  J0 


rr  [*■<*( 


(if)  (2wxe  +  vx)£'j  dx  tie 


2Lf  f2ff  fw  Tue  /w0+v\z  ?  ,  ue  /  ue\1 

J0  \^Li^+ (~R~)  +  v(w*+v*)'|1’v)ir(v*+R‘)j 

[.W  +  Vni  wn  +  V  w  Un  *  ^ 

2vw*(-r ~)  +  (i;v)(-r- )(v*+t)  +  2u*>] 


2(  1  -v£)  *'0  *'0 
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It  will  be  assumed  that  the  buckling  displacements  can  be  represented  by 
Equation  [10]: 


«  " *  I  /V-yj  A. 

u  =  cos  n0  2^  “m  coa|-~_J 


8in  ”9  Vm  3in  f  -T7~  : 

\  i-V  / 

m=  1 


“  5  v 

0  v*  •  i  m  \ 

^  ™  \  R  ) 

'in-- 1 


38 


Before  proceeding  with  the  integration,  a  few  assumptions  can  be  made 
which  simplify  the  results, considerably .  These  are  as  follows: 


»  v 


A2ir  ^  J*2  7T 

Wq  d.8  »  |  v.q  d$  ,  j  R  U0  v.,f  cl 8 
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wjde  »  I  v^de  ,  j  — ■  de 

Jo  fo  Jo  R 
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v2  ae 


The  ba#is  for  these  conditions  is  the  assumption  that  the  relative  orders 
of  magnitude  of  vm,  and  wm  can  be  estimated  from  the  results  of 

the  stress  function  approximation  (Equations  [12]).  The  assumptions  in 

through  [27],  and  are  based  bn  the  fact  that  unity  is  negligible  compared 
with  n2  and  A^  for  short  shells.  If  the1  small  quantities  in  Equation  [49] 
are  neglected,  Equation  [43]  is  reduced  to 
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When  the  series  expressions.  Equations  [6]  and  [10],  are  substituted 
into  Equation  [50]  and  the  integrations  are  performed,  the  reault  is 
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where,  for  convenience,  the  equation  has  been  divided  by - = —  . 

(1~v4)R 

The  boundary  conditions  (Equation  [47])  on  w  and  v  are  automatically 
satisfied  by  Equations  [10];  the  conditions  on  wx  and  u0  require 
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[52] 
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so 
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Thus  summations  of  the  odd  and  the  even  terms  must  vanish  separately. 

The  problem  now  is  to  arrive  at  a  stationary  value  of  the  change  ip 
total  potential,  Equation  [SI  ],  while  at  the  same  time  satisfying  the  four 
conditions  of  constraint,  Equation  [53].  This  can  be  done  through  the 
use  of  the  Lagrange  multiplier  method  which  requires  that  the  expression 
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be  stationary,  where  the  cr's  are  the  unknown  multipliers.  The  require¬ 
ment  gives  rise  to  the  equations 


8U1  _  8U^  =  MT 

8uj  ”  9vj  “  9Wj 


[55] 


which  can  be  written  explicitly  as  follows:' 


A.  U.  +  Bj  v..  +  C.  w.  +  5"  c, ..  w.  +  —  [l  -  (-l)"1]  +  [l  +  ( - 1  )’h  K  0 

J  .1  .5  J  J  J  J1  1  2  1  1  2  1 

i*  1 
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ff, 
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The  quantities  designated  by  Aj,  Bj;  Cj .  etc.  ,  should  not  be  confused 
with  other  quantities  appearing  earlier  in  mis  report  which  bear"  the  same 
notation. 
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SjUj  +  Dj  Vj  +  Ejwj  +  2,e)iW>  =  ° 
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It  will  be  seen  that  once  again  the  coupling  conditions  on  the  odd  and  the 
even  modes  are  in  effect  because  of  the  coefficients  multiplying  the  ar's 
and  the  definition  of  am.  Thus  Equation  [56]  represents  two  systems  of 
equations  corresponding  to  all  odd  and  all  even  values  of  j  and  i.  As  be¬ 
fore,  only  the  odd-valued  system  is  of  practical  interest.  If  j  and  i  are 
given  an  upper  limit  J,  the  complete  set  bf  equations  for  this  system  is 
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the  last  two  equations  representing  the  two  constraints.  There  are  thus 

2  +  Aliiii  equations  and,  including  and  a^,  2  +  unknowns. 

The  system  can  be  solved  in  the  same  manner  as  Equation  [32],  How¬ 
ever,  it  is  highly  complicated;  moreover,  because  c-.  and  eJ5  are  not 

J*  J*- 

generally  invariant  when  i  and  j  are  interchanged,  the  matrix  of  the  co¬ 
efficients  is  not  symmetric. 

This  situation  can  be  greatly  simplified  if,  as  in  the  case  of  Equa¬ 
tion  [34] ,  we  neglect  all  c--,  e--7  and  f--  for  which  j  and  i  are  unequal. 
Further  simplifications  can  also  be  made  by  neglecting  the  small  quantity 
4>(1  -  v  )  appearing  in.  the  coefficients  B •  and  Cj,  and  -1  in  the  second 
term  of  Fj.  Equations  [57]  then  reduce  to 
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The  equations  can  be  solved  in  a  manner  similar  to  that  followed 
in  the  solution  of  Equations  [  35  j .  From  the  first  three  equations,  vj 
ia  eliminated  and  the  expressions 
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are  obtained.  When  these  are  substituted  into  the  last  two  of  Equations 
[58]  and  or j  and  a ^  are  eliminated,  the  resulting  buckling  equation  is 
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By  neglecting  the  small  terms  involving  i(t  ,  the  complexity  of  the  expres¬ 
sions  is  greatly  reduced,  and  after  some  rearrangement,  the  equation 

takes  the  form 
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The  solution  of  Equation  [61]  can  be  accomplished  graphically  with  the 
game  procedure  described  for  Equation  [39]>  Once  again  the  lowest  root 
should  occur  between  the  first  two  asymptotes  where  Oj  and  Q-j  vanish. 
The  accuracy  of  the  result  will  depend,  as  before,  on  the  size  of  J. 


NUMERICAL  EXAMPLE 


"-Tfeir "tsswraplif tirb* used  fbr numerical  study  is  the  test  of  the  ma¬ 


chined  cylinder  BR-4B  having  external  rectangular  stiffeners  and  which 
is  considered  in  the  section  on  experimental  work.  Its  dimensions  in 
inches  are  as  follows: 


R  ~  8.048 
h  a  Q. 081 
Lf  -  4.266 
d  (frame  depth)  *  0.570 
b  (frame  width)  =  0.  138 
p  =6.80 

The  value  of  p  (6.80)  is  considerably  larger  than  that  usually  encountered 
in  pressure-vessel  design. 

Rings  of  Finite  Rigidity 

As  a  preliminary  to  the  calculations,  Table  1  shows  the  first  11  co¬ 
efficients  (a  )  of  the  initial  deflections  (w  and  u)  as  calculated  by  Equa¬ 
tion  [8];  m  is  written  with  the  absolute  value  sign  since  both  positive  and 
negative  values  are  permitted.  Despite  the  relatively  large  value  of  p, 
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Table  2 


for  Cylinder  BR-4B 


i 

!  i  * 

i  lm! 

S 

s 

I 

* 

f 

am  j 

j  Equation  [8]  j 

1  0 

j  1.00000 

1 

|  0. S4590 

2 

0.25538 

3 

0.06345 

4 

0. 02099 

5 

0.00870 

1 

/ 

6 

0.00422 

7 

0, 00228  J 

8 

0.00134 

9 

0.  00084 

10 

0.00055 

11 

0.00037 

Comparison  of  Approximate  and 
Matrix  Solutions  for  Cylinder 
42  with  5^ ini t s  H.’’‘r'rTO 
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i  - — — — — 

\  Buckling  Pressure,  osi 

\ - - 

i 

Matrix 
Equation  [33] 

A.pprcxim:  e  j 
Equation  |3Qj 

1 

693.  I* 

AC31*5*  ! 

~ '  i 

3 

657.1 

663 

5 

646.3 

654 

n 

« 

641.4 

648 

Q 

f 

638.7 

645 

11 

•636.9 

642 

13 

635.7 

Calculations 

i 

A84  8 

performed  on 

IBM-  7  090 

17 

634.  1 

computer 

19 

633.5 

^Slide  rule 

21 

633.  1 

calculations 

Pressures  for  n=  11,  , 

E  =  28. 9  x  10b  psi 


— 

the  convergence  of  23  am  *B  rapid.  ror  purposes  of  calcula¬ 

te  -a 

tion,  only  the  first  four  or  five  values  are  significant.  As  previously  indi 
cafedi  the  convergence  wili’be  more  rapid  for  smaller  values  of  p. 

It  is  of  some  interest  to  compare  results  of  the  matrix  solution. 
Equation  [33],  and  the  simplified  solution.  Equation  [39].  The  critical 
buckling  mode  for  both  equations  was  found  to  be  11.*  Table  2  shows 
the  buckling  pressures  according  to  these  two  solutions  for  successive 
values  of  J.  Since  Equation  [33]  was  solved  using  an  IBM-7090  com¬ 
puter,  results  were  readily  obtainable  for  large  values  of  J  {an  arbi¬ 
trary  limit  of  21  was  imposed).  Solutions  of  Equation  [39]  were  obtained 


Before  the  calculations  were  begun,  a  preliminary  estimate  of  the  crit 
cal  value  of  n  was  obtained  from  Figure  8  for  the  case  of  simple  support. 
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a: 


by  slide  rule  up  to  J  =  11;  it  did 
not  appear  worthwhile  to  continue 

beyond  this. 

The  information  in  Table  2  is  **° 

also  shown  graphically  in  Figure  10.  |.  i7° 

It  is  evident  that  the  convergence  of  f  ** 
both  solutions  is  fairly  rapid  so  that  f  **° 
little  accuracy  is  gained  in  proceed-  j  **° 
ing  beyond,  say,  J  =  7.  Such  a  con-  ^ 
elusion,  however,  cannot  be  stated  *° 

generally.  If  the  frames  had  been 
larger,  the  solutions  would  probably 
have  converged  more  slowly.  For 

...  ,  ..  .  ..  .  Figure  10  -  Variation  in  Approximate 

this  example  it  also  appears  that  and  Matrix  Solutions  with  J 

Equation  [39]  is  a  reasonably  ac-  ^  .'for  Cylinder  BR-4B 

curate  approximation  to  Equation  ' 

[33].*  Furthermore,  since  the  accuracy  of  the  approximation  depends 

CD 

only  on  the  rapidity  of  convergence  of  £  24^ ,  one  should  expect  closer 

m=~® 

agreement  for  smaller  values  of  p. 

Equation  [33]  also  permits  the  determination  of  the  relative  ampli¬ 
tudes  wj  of  the  buckling  modes.  Table  3  lists  the  ratios  Wj/w^  for  n  =  11 
and  J  =  21.  It  will  be  noted  that  the  components  are  of  successively  de¬ 
creasing  magnitude  and  the  first  is  much  larger  and  cf  opposite  eign  from 
all  the  rest.  The  rotation  of  a  frame  wi^i/r  produced  by  the  first  com¬ 
ponent  is  thus  partially  reduced  by  the  sum  of  the  succeeding  components. 

Simple  Support  Conditions 

Results  of  calculations  for  Cylinder  BR-4B  (j3=  6.80}  with  simple 
supports  replacing  the  ringe  are  given  in  Table  4.  The  pressure  and  mode 
obtained  from  the  minimized  curves  cf  Figures  8  and  9,  representing  Equa 
tions  [43]  and  [44],  compare  well  with  the  results  for  the  unminimized 

*Since  Equation  [39]  does  not  contain  the  approximation  in  the  linear 
case,  the  results  of  both  equations  should  be  identical  for  J  =  i. 


Equation  [42] „  Corresponding  results  for  the  Von  Mises  solution,  also 
obtained  from  Figures  8  and  9,  show  a  slightly  higher  pressure  and  dif¬ 
ferent  mode.  This  difference  reflects  the  influence  of  the  boundary  con¬ 
ditions  on  the  initial  deformations.  Although  for  this  case,  the  inclusion 
of  this  effect  gives  rise  to  a  lower  pressure.  Figure  9  shows  that  the 
pressures  will  be  higher  for  fi  <  4.55  , 

Rings  of  Infinite  Rigidity 

Using  Equation  [hi  J,  it  should  be  possible  to  approach  aa  closely 
as  desired  the  exact  solution  to  the  problem  of  fixed  support,  consistent 
with  the  approximations  made  in  the  derivation  of  the  equation.  It  may 
be  surprising,  then,  to  consider  the  results  of  calculations  for  Cylinder 
BR-4B  given  in  Table  5.  Equation  [61],  carried  out  to  the  eleventh  degree 


'  Table  3  • 

Relative  Amplitude*  of  Buckling 
Modes  for  Cylinder  BR-4B 
with  Finite  Rings 

i — ' — \  - - — — : - — i 
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Ttble4 

Buckling  Pressures  for  Cylinder 
BR-4B  with  Simple  Supports 


I  "V 


t — m 


i 

Equation  [33] 
nail  J  =21 

i 

i. 00000 

3 

-  0.06654 

5 

-  0.01962 

7 

-  0.00773 

9 

-  0.00378 

11 

-  0.00211 

13 

-  0.00129 

15 

-  0. 00085 

17 

-  0.00058 

19 

-  0. 00042 

21 

-  0.00031 

Equation  [42] 

550 

11 

Minimised  Form 
(Figures  8  and  9) 

551 

11 

Minimised  Von  Mines 
Solution!39 
(Figures  8  and  9) 

580 

12 

Shell  length  assumed  to  be 
ccnler-to-center  frame  spacing 
(Lf-~4.266  in.) 

Table  5 


Buckling  Pressures  for  Cylinder  BR-4B 
with  Infinitely  Rigid  Rings 


Effect,  of  Boundary 
Conditions  onj 

Initial  Deflections 

Theory 

Pc 

psi 

n 

Included 

[*"’  *  ,Ti7ip] 

Equation  [6l] 

(J  =  21) 

811 

H 

Neglected 

Equation  [61  ] 

(J  «  21) 

* 

810 

11 

*  °) 

Nash38 

804 

12 

In  all  cases,  shell  length  was  taken  to  be  the  eenter- 

to-center  frame  spacing  ( L*  ~  4  266  in. ) 


(J  =  21),  gave  a  buckling  pressure  of  811  psi  for  n=  11,  while  the  far  sim¬ 
pler  equation*  of  Nash9*  gave  804  psi  for  ns  12.  This  latter  equation  is 
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agreement  with  other  more  complicated  solutions  also  obtained  by  Nash.16 
The  fundamental  characteristic  of  all  of  these  solutions  is  that  they  are 
based  on  the  assumed  budding  configuration 


W  »  W„ 
G 


1 


[62] 


Eh  f 

' R  X 


12R2(l-v2} 


4  u. 


8X]  + 


.  ..2  2 

+  4X  j,  n 


8X1  ] 
(4\2  +  n2j2  / 


-  (n2  -  1)  +  X2 

2  1 


In  all  calculations,  the  distance  between  supports  wag  taken  to  be  the 
center -to -center  frame  spacing  ( Lf  =  4.266  in.).  Frequently,  when  for¬ 
mulae  of  this  type  arc  applied  to  ring  -  stiffened  cylinders,  the  inner  or 
unsupported  distance  between  frames  is  \saed. 
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This  function  will  not,  in  general,  satisfy  the  differential  equations  of 
equilibrium.  Rather  it  has  been  chosen  somewhat  arbitrarily  because  it 
satisfies  the  boundary  conditions  and  is  conveniently  used  with  the  energy 
method.  When  applied  to  ring  stiffened  cylinders,  the  nolutions  based  on 
this  function  have  often  predicted  pressures  much  higher  than  those  ob- 
served  experimentally.  Consequently,  cm®  natural  conclusion  is  that  the 
use  of  a  different  function  could  lead  to  lower  pressures.  While  the  re¬ 
sults  in  Table  5  are  too  limited  to  warrant  a  general  conclusion,  they  in¬ 
dicate  that  the  configuration  of  Equation  [62]  may  in  fact  be  ne&sly  correct 
and  that  the  difficulty  actually  lies  in  the  assumption  that  full  fixity  can  be 
provided  by  rings  of  practical  sise. 

It  is  also  worth  noting  that  the  convergence  of  Equation  [61]  was 

much  slower  than  for  the  other  two  cases  previously  discussed,  This  was 

J 


caused  by  the  slow  convergence 


of  r 

ij=i 


h 


Q 


j 


which  necessitated  carrying 


the  calculations  out  to  .T  =  2 1 :  From,  this  it  may  be  concluded  that  the  sine 
series  (Equation  [10])  is  not  well  suited  to  the  boundary  conditions  for 

for  example,  to  Equation  [62],  would  probably  result  in  more  rapid  con¬ 
vergence. 

Table  5  also  shows  that  for  the  seometrv  in  ni»*Hn»i  ♦**  - 

......  f  vwuuiuuiiii  uu  uie  imuai  oenections  had  a  negligible  effect 

on  the  collapse  pressure.  When  this  influence  is  neglected  (am  =  0),  the 
initial  deflections  are  those  for  an  infinitely  long  cylinder,  the  condition 
assumed  by  Nash,  and  the  buckling  pressure  is  810  psi.  That  the  pres¬ 
sure  (804  psi)  given  by  tine  Nash  aquation  is  lower  than  this  can  probably 
be  explained  by  the  fact  that  the  buckling  configuration  assumed  by  Nash 
does  not  completely  satisfy  the  boundary  conditions  of  Equation  [47] .  Be¬ 
cause  of  certain  approximations  arising  from  the  use  of  a  stress  function 
of  the  Donnell  type,  the  circumferential  buckling  strain  does  not  van¬ 
ish  at  the  boundaries  * 


r  ■'t  ’’X 


^Following  the  procedure  indicated  in  Appendix  E,  Nash  eliminates  u 
and  v  through  the  use  of  a  stress  function.  Because  w  for  this  case  is 
of  the  form  of  Equation  [62],  the  resulting  stress  function  is 
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In  view  of  these  results,  it  might  be  of  interest  to  compare  the  buck¬ 
ling  configuration  employed  in  the  present  analysis  with  that  represented 
by  Equation  [62]  for  this  example  on  a  quantitative  basis.  By  suitable  re¬ 
arrangement  of  Equations  [58]  through  [oil,  an  equation  expressing  the 
relative  magnitude  of  the  corresponding  buckling  coefficients  Wj  is  ob¬ 
tained: 


w- 


Wi 


[63] 


where 


j»k  *  1,  3,  5 . J 


,^,,|{l^S|^KCOft|im*ed  (torn preceding  ... 

ER\| 

F»-  t  ■  "  w  cos  n0  cos 
n+X“ 

in  the  notation  of  this  report,  and  the  . . . 

1  /82F  v82F  \ 

°9 


does  not  vanish  at  x  =  0  and  x  *  .  However,  if  Equations  [122]  of 

Appendix  E  are  modified  so  that 


<*x  » 

(JQ  - 

^xe  = 


82F 


+  f(0) 


R2892 

i  g(e) 


82F 


8x4 

82F 
R  8x  88 


it  is  then  possible  to  satisfy  (he  condition  on  by  proper  adjustment  of 
the  functions  f(9)  and  g(0).  The  resulting  buckling  pressure  is  slightly 
higher  than  that  given  by  Nash. 
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and1  all  quantities  retain  the  definitions  of  Equation  [  61] .  An  equivalent 
of  Equation  [  63]  can  be  obtained  for  the  buckling  configuration  of  Nash 
by  expanding  the  function1^ 


x  \ 


w  =  w„  I  1  -  cos  -- —  } 
°  \  R  / 


0  S  X  £  L  t 


L< 


W  =  -  WQ  j  1  -  COS 


£  x  £  2L>£ 


[64] 


in  a  Fourier  scries.  The  result  1b 


16wq 

*j(jZ-4> 


J'1,3,5 


[65] 

J 


In  order  to  compare  these  results,  calculations  were  carried  out 
withthe  aid  of  Equation  [  63]  for  flie  humer ical  example  where 
J  «  21,  and  pc  -  810  pai.  The  results  appear  in  Table  6  in  the  form  of 
Wj/wj,  together  with  the  corresponding  values  from  Equation  [65].  Ratios 
bevond  I  s  11  were  not  calculated.  Tt  will  h*  naan  that  thum  {*  ron§i dur¬ 
able  difference  in  the  individual  coefficients.  Apparently  their  total  effect 
is  such  that  the  difference  in  buckling  pressures  is  insignificant. 

SUMMARY  AND  CONCLUSIONS 

1.  An  analysis  for  the  elastic  buckling  of  ring-supported  cylindrical 
shells  has  been  developed  which  considers  the  influence  of  the  rings  on 
deformations  before  and  during  buckling, 

2.  The  bending  stresses  associated  with  the  prebuckling  deforma¬ 
tions  have  a  negligible  effect,  on  the  theoretical  buckling  strength. 


^Tha  function  in  Equation  [64]  is  periodic  in  the  interval  0  <  x  <  2Lf 
while  Equation  [62]  repeats  itself  in  each  bay.  However,  their  corres¬ 
ponding  buckling  pressures  are  identical. 


53 


Table  6 


Comparison  of  Buckling  Configurations  for  Cylinder 
BR-4B  with  Infinitely  Rigid  Rings 


j 

- - -  ■  . . — . ■-  . . 1  g 

Wj/W, 

Equation  [63]* 

Equation  [65]  ■ 

(Nash35) 

1 

1. 0000 

1 . 0000 

3 

-0, 1492 

-0.2000 

5 

-0. 0387 

-0.0286 

mm 

-0. 0173 

-0.0095 

B 

-0.0085 

-0.0045 

m 

-0. 0047 

•U, wvfcj 

tm  =  0 ;  J  =  2 1  ;  n  =  11  ;  pc 

s  810  psi. 

3.  The  buckling  equation*  for  a  shell  with  finite  rings  represented  by 
the  stability  determinant,  Equation  [33],  can  be  approximated  with  good 
accuracy  by  the  single  Equation  [39]. 

4.  A  special  result  of  the  general  analysis  is  the  buckling  equation 
for  a  shell  with  simple  supports.  This  differs  from  the  buckling  equation 
of  Von  Mises  since  it  accounts  for  the  effect  of  the  boundary  conditions  on 
the  prebuckling  deformations . 

5.  Anothe”  special  case  is  that  for  a  shell  with  fully  fixed  edges.  It 
appears  possible  from  numerical  results  that  the  use  of  tine  simplified 

i% 

form  1  -  cos  ~ —  to  represent  the  longitudinal  budding  profile  may  be  a 
Lf 

reasonable  approximation  for  this  case. 

6.  The  energy' approach  has  been  used  to  obtain  the  trigonometric 
series  equivalent  of  the  solution  of  Pulos  and  Salerno  for  the  initial  axi~ 
symmetric  deformations  of  a  ring -stiffened  cylinder.  This  form  has  cer¬ 
tain  mathematical  advantages,  particularly  when  employed  in  the  solution 
of  buckling  problems. 
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PART  n  -  EXPERIMENT 


EARLIER  TESTS 


Following  the  early  work  of  Windenburg,29  later  experimental  studies 
at  the  Model  Basin  were  directed  toward  the  evaluation  of  classical  buck¬ 
ling  theory  for  snorj,  shells.  Tests  of  ring- stiffened  cylinders,  designated 
BR-1'°  and  BR-531,  which  were  fabricated  from  steel  plate  by  welding  and 
not  stress  relieving,  resulted  in  collapse  pressures  well  below  theoreti¬ 
cal  expectations.  It  was  suspected  that  these  discrepancies  could  be  at¬ 
tributed  to  weakening  effects  of  imperfect  circularity  and  residual  rolling 
and  welding  stresses.  To  investigate  this  suspicion  further,  two  additional, 
cylinders,  BR-4  and  BR-4A  ,  were  tested.  These  were  identical  in 
geometry  and  in  material  properties,  but  BR-4  was  rolled  and  welded 
while  BR-4A  was  machined  from  a  stress -relieved  thick  cylinder.  The 

f!  h  4  ft  T2.  r\fVi  ttrA« m  tt  ri  f  1  /xKao 
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appearing  between  stiffeners.  However,  BR-4  collapsed  at  390  psi  where¬ 
as  BR-4A  collapsed  at  550  psi. 

With  the  weakening  effects  of  imperfections  so  clearly  demonstrated, 
it  was  evident  that  further  attempts  to  obtain  elastic  buckling  data  with  im¬ 
perfect,  fabricated  cylinders  would  be  fruitless.  A  better  alternative  ap¬ 
peared  to  be  the  study  of  machined  and  initially  Stress-free  cylinders. 

This  approach  had  already  proven  successful  in  experimental  studies  of 
general  instability  Accordingly,  a  third  test  cylinder,  designated 

BR-4B,  was  manufactured.  Geometrically,  this  was  a  smaller  scale  ver¬ 
sion  of  BR-4  and  BR-4A,  but  it  was  machined  from  tubing  of  higher  yield 
strength.'  With  this  combination  of  properties,  it  was  expected  that  col¬ 
lapse  would  be  initiated  by  elastic  buckling. 

DESCRIPTION  OF  TEST  CYLINDER 

BR-4B,  whose  dimensions  were  given  as  a  numerical  example  ear¬ 
lier  in  this  report,  was  a  b/10-scale  model  of  BP.-4A.  As  shown  in  Fig¬ 
ure  11,  it  consisted  of  four  central  full-length  bays  and  a  short  bay  !at 
each  end  terminated  by  a  heavy  bulkhead  ring.  One  ring  had  evenly  spaced 
holes  to  accept  bolts  for  a  flat  closure  head.  The  four -bay  arrangement 
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was  intended  to  provide  a  central 
test  section  free  of  end  effects.  The 
cylinder  was  machined  from  a  thick, 
forged- steel  tube  which  had  been  part 
of  a  gun  barrel  liner.  Compression 
tests  cf  specimens  taken  from  the 
tube  prior  to  machining  indicated  a 
yield  strength  of  82,500  psi  {based 
on  an  offset  strain  of  0.002}  com¬ 
pared  with  50,000  psi  for  BR-4  and 
ER-4A,  The  elastic  limit  was  found 
to  be  74,400  psi.  Since  the  elastic 
buckling  strength  of  the  cylinder 
would  be  directly  proportional  to  the 
modulus  of  elasticity  of  the  material, 

w  A  ...«  _  _  A  i.t_  .  A.  i.L  i  •  V.  . 
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determined  accurately.  For  (his  pur¬ 
pose,  a  ring  specimen  was  taken  from  one  of  the  bulkhead  rings  at  the  con¬ 
clusion  of  the  tests.  Measurements  of  the  deflections  of  the  ring  under 
diametral  point  loading  were  used  to  establish  a  modulus  of  28.9  x  1()6  psi. 
This  procedure  is  described  in  Appendix  F. 

PRELIMINARY  TEST 

One  of  the  objectives  of  the  investigation  was  to  determine  whether 
any  region  in  the  area  of  failure  had  been  stressed  beyond  the  elastic  limit 
prior  to  collapse.  Ideally  this  would  indicate  extensive  strain-gage  instru¬ 
mentation  of  this  critical  area,  but  to  locate  such  an  area  prior  to  testing 
if!  an  impossibility.  On  the  other  hand,  it  was  not  practical  to  .provide  each 
bay  with  such  extensive  strain-gage  instrumentation  that  adequate  coverage 
of  any  area  would  be  insured.  Consequently,  it  was  decided  that  a  prelim¬ 
inary  test  should  be  carried  out  in  an  effort  to  determine  in  which  bay  buck¬ 
ling  would  initiate.  If  this  could  be  accomplished,  the  critical  bay  would 
then  be  liberally  instrumented  and  a  second  test  carried  out  in  which  the 
cylinder  would  be  taken  to  collapse.  An  additional  possibility  vv<*s  that 
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Figure  11  -  Test  Cylinder  BR-4B 
All  dimensions  in  inches 
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nonlinear  elastic  strains  would  appear  during  the  preliminary  tests,  and 
would  be  of  such  a  magnitude  that  the  buckling  pressure  could  be  predeter¬ 
mined  by  means  of  the  Southwell  method.  This  procedure  had  been  used 
with  considerable  success  in  tests  of  cylinders  collapsing  by  general  in¬ 
stability,38  but  as  yet  had  not  been  tried  in  the  case  of  shell  instability. 
Since  nonlinear  strains  were  not  observed  in  the  test  of  BR-4A,  it  was  not 
expected  that  they  would  appear  at  pressures  below  540  psir  the  highest 
pressure  at  which  strains  were  measured  dttring  that  test. 

The  exterior  of  the  cylinder  was  instrumented  with  electrical  resis¬ 
tance  strain  gages  located  circumferentially  at  intervals  of  10  deg  in  the 
middle  of  each  of  the  four  typical  bays  and  extending  completely  around 
the  circumference.  Temperature  compensation  was  provided  by  "dummy” 
gages.  The  test  chamber  was  a  ?,500-psi  pressure  tank,  20  in.  in  diame¬ 
ter,  and  oil  was  used  as  the  pressurizing  fluid.  Prior  to  the  actual  test, 
the  cylinder  was  immersed  in  the  tank  in  a  free-flooding  condition.  The 
pressure  was  then  raised  in  increments  to  500  psi  and  strains  measured 
to  detect  gages  that  were  undesirably  "pressure-sensitive.  "  Those  which 
exhibited  a  sensitivity  of  5  pin. /in.  or  more  for  100  psi  were  considered 
unsatisfactory.  These  gages  were  then  checked  by  observing  the  strain 
induced  when  local  pressure  was  applied  to  the  gage.  It  was  found  that 
this  procedure  was  an  adequate  substitute  for  the  pressurizing  method, 
and  it  was  used  subsequently  as  gages  were  successively  replaced  and 
checked  until  all  were  satisfactory. 

The  cylinder  was  then  placed  in  the  tank  with  one  and  closed  by  a 
fla*,  circular  plate  and  the  other  sealed  against  the  tank  top,  which  had 
an  opening  to  permit  access  to  the  interior  of  the  cylinder.  While  strains 
were  read  with  automatic  strain  recorders,  the  pressure  was  applied  in 
small  increments  up  to  570  psi.  Ail  strains  were  still  linear  at  this  pres¬ 
sure,  and  it  appeared  that  proceeding  to  a  higher  pressure  would  only  in¬ 
crease  the  risk  that  the  cylinder  would  collapse  before  more  extensive  in¬ 
strumentation  could  be  installed.  The  pressure  was,  therefore,  reduced 
in  increments  back  to  zero. 

It  was  apparent  from  the  strain  plots  that  no  yielding  had  taken  place. 
The  largest  apparent  permanent  set  for  all  operating  gages  was  45  pin.  / in. 
which  was  considered  insignificant  in  a  total  strain  of  1500  to  1800  pin. /in. 


Table  7 


Circumferential  Strain  Sensitivities  (Exterior) 
at  Midbay,  Preliminary  Test 


Bay 

- 

Compressive  Strain  Sensitivity  in  j 

micro-inches  per  inch  per  psi  j 

Maximum 

Minimum 

Average 
of  All 
Gages 

l 

Z.  91 

2.52 

2 

2.96 

2.40 

3 

3.  15 

2.72 

4 

3.  12 

2.79 

Theory  of  Reference  33:  Sensitivity  = 

3.03 

and' probably  due  to  the  overall  effect  of  pressure  sensitivity  and  drift  in 
the  measuring  system.  The  strain  sensitivities  in  micro-inches  per  inch 
per  pound  per  square  inch  observed  during  the  tests  are  summarised  in 
Table  7.  The  average  sensitivity  was  highest  (2.97)  in  Bay  3  and  next 
highest  (2.92)  in  Bay  4.  The  maximum  individual  sensitivity  was  observed 
at  90  deg  in  Bay  3.  On  the  basis  of  these  measurements,  it  appeared  that 
Bay  3  woixld  be  critical,  but  that  Bay  4  would  also  deserve  close  attention. 
Final  instrumentation  of  the  cylinder,  carried  out  on  the  basis  of  these 
observations,  is  described  in  the  next  section. 

FINAL  TEST 

Instrumentation  and  Test  Procedure 


Gage  locations  for  the  final  tost  are  shown  in  Figure  12.  This  time, 
interior  a3  well  as  exterior  gages  wore  installed,  the  total  number  being 
limited  to  the  capacity  of  Lhe  available  automatic  recorders.  Most  u l  the 
gages  were  concentrated  in  Bays  3  and  4  with  major  emphasis  on  Bay  3. 
The  pairs  of  circumferential  and  longitudinal  gages  were  sc  arranged  that 
stresses  could  be  measured  at  the  two  locations  { outside  at  midbay  and 
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Figure  12a  -  Bays  1  and  2 
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Figure  12b  -  Bays  3  and  4 
Figure  iZ  -  Gage  Locations  for  Final  Test 
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inside  at  a  frame)  where  they  are 
normally  highest,  The  midfcay  ar¬ 
rays  of  circumferential  gages  were 
located  in  the  four  typical  bays  for 
the  purpose  oi  detecting  iobe  pat¬ 
terns  should  they  appear  prior  to 
collapse.  These  measurements 
would  be  utilized  in  determining 
the  elastic  buckling  pressure  by 
means  of  the  Southwell  method.5* 
The  exterior  circumferential  gages 
located  in  Bay  3  along  the  90-deg 
generator  were  intended  to  give  an 
indication  of  the  longitudinal  pro¬ 
file  of  a  lobe  in  the  event  one  ap¬ 
peared.  Since  the  highest  strain 
in  the  preliminary  test  was  meas¬ 
ured  at  90  deg  In  Bay  3,  this  appeared  to  be  a  likely  location  for  the  devel¬ 
opment  of  a  lobe  trough.  Gages  were  also  mounted  on  the  aides  of  Frame 
4  to  provide  a  meaaure  of  any  twisting  and  bending  that  might  take  place 
prior  to  collapse.  All  exterior  gages  were  examined  for  pressure  sensi¬ 
tivity  and  replaced  in  &ec  o i danc s  with  the  procedure  followed  for  the  pre¬ 
liminary  test. 

Two  pressure  runs  were  made  during  the  final  test,  as  indicated  in 
Table  8.  On  the  first  run,  a  maximum  pressure  of  610  psi  was  attained. 
During  the  second  run,  collapse  occurred  suddenly  at  633  psi,  the  last 
strain  readings  having  been  made  at  630  psi,  and  the  pressure  immediate¬ 
ly  dropped  to  a  residual  value  of  377  psi. 

Results  (iiiu  Discussion 


Table  8 


Loading  Schedule  for  Final  Test 


Pressures  (psi)  at  which 

Strains  Were  Measured 

First  Run 
■■■  — * — 

Second  Run 

0 

570 

0 

610 

50 

580 

100 

620 

100 

590 

200 

630 

200 

600 

300 

Failure 

300 

610 

400 

(633) 

400 

570 

500 

377 

500 

500 

540 

(Residual) 

520 

300 

560 

540 

100 

580 

0 

560 

0 

600 

Visible  damage  to  the  cylinder,  as  shown  in  Fig\ire  13,  was  con¬ 
fined  to  Bays  3  and  4  where  a  typical  lobar  pattern  developed  in  the  area  of 
heaviest  instrumentation.  While  the  lobes  did  not  completely  encircle 
the  circumference,  a  close  examination  indicated  that  a  complete  pattern 
would  have  contained  !  I  lobes. 
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Figure  1 3a  -  Exterior  View  Figure  i3b  -  Interior  View 

Figure  13  -  Views  of  Collapsed  Cylindor  BR-4B 

During  the  first  run,  the  strains  were  generally  linear  with  pressure 
up  to  about  590  psi.  Small  deviations  from  linearity  appeared  at  600  psi, 
and  became  more  noticeable  when  610  psi  was  reached.  On  the  unloading 
portion  of  the  run,  the  strain-pressure  plots  obtained  during  loading  were 
in  moBt  cases  retraced.  At  the  conclusion  of  the  run,  the  maximum  zero 
shift  was  approximately  50  pin.  / in.  There  appeared  to  be  no  correlation 
in  general  between  zero  shifts  and  departures  from  linearity  observed 
during  the  pressure  run.  As  in  the  preliminary  run,  the  shifts  can  be 
attributed  to  error  in  the  measuring  system  rather  than  to  a  significant 
yielding  of  the  material. 

On  the  second  run,  the  deviations  from  linearity  again  appeared  at 
600  psi  and  grew  with  pressure  until,  at  630  psi,  a  pronounced  lobar  pat¬ 
tern  was  evident.  Such  patterns  have  often  been  observed  in  stiffened  cyl¬ 
inders  prior  to  collapse  by  general  instability.38  They  are  produced  by 
initial  imperfections  which,  even  though  minute  in  machined  cylinders, 
have  a  decided  effect  near-  the  point  of  budding,  Figure  14  shows  some 
examples  of  the  strain  pressure  plots  for  both  runs.  The  circumferential 
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Figure  14  -  Typical  Plots  of  Circumferential  Strains 
at  Midbay  versus  Pressure  (Final  Test) 

lobar  strain  patterns  in  Bays  3  and  4  at  630  psi  are  shown  in  Figure  15, 
For  contrast,  the  strains  at  500  psi  are  also  shown  It  will  be  seen  that 
the  patterns  in  the  two  bays  are  staggered,  just  as  are  the  final  deforma¬ 
tions  in  Figure  13,  Although  the  strain  patterns  are  not  complete  because 
of  a  few  gage  failures,  close  inspection  of  both  patterns  indicates  the  ex¬ 
istence  of  13  lobes.  It  can  also  be  seen  that,  as  anticipated,  a  lobe  trough 
developed  at  or  near  the  90-deg  generator  of  Bay  3. 

From  strain  data  of  the  type  shown  in  Figure  14,  it.  was  possible  to 
obtain  Southwell  plots  of  fairly  good  quality.  Plots  for  the  gages  of  Figure 
3  4  appear  in  Figure  16.  With  this  method,  the  nonlinear  component  («*) 
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Figure  15  -  Circumferential  Strain  Patterns 
(Run  Z,  Final  Test) 


of  the  total  strain  is  determined  and  plotted  against  the  ratio  -jj—  .  T  The 
elastic  buckling  pressure  is  then  obtained  from  the  slope  of  the  resulting 
straight  line  plot.  This  was  done  for  a  number  of  midbay  circumferential 
gages  located  externally  in  Bays  3  and  4,  The  gages  were  selected  on  the 
basis  of  data  quality  from  Jiose  showing  the  largest  deviations  from  line¬ 
arity.  Table  9  gives  a  summary  of  these  results.  It  can  be  seen  that  the 
pressures  are  all  in  good  agreement.  Furthermore,  the  average  buckling 
pressure  of  637  psi  in  very  close  to  the  experimental  collapse  pressure 


t Since  c*  is  positive,  or  negative  depending  on  the  gage  orientation, 

It  t3?rj 

and  ”|  were  used  in  Figure  16  so  that  all  points  would  lie  in  the 
first  quadrant. 
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Table  9 


Figure  16  -  Typical  Southwell  Plots 
for  Exterior  Circumferential 
Strains  at  Midbay 
(Run  2,  Final  Teat) 


Southwell  Results  from 
Circumferential  Strains 
at  Midbay 

(Run  2,  Final  Test) 


Bay 

Angular 

Location 

B 

3 

H 

639 

636 

634 

180° 

643 

4 

40* 

638 

60® 

636 

90® 

638 

no® 

634 

2L  r'znr -f  * 

. .. . 

639 . 

180® 

_ 

636 

Average  pc  =  637 

Gages  selected  on  basis 
of  data  quality  from  those 
showing  largest  nonlinear 
strains. 

of  633  psi.*  However,  it  should  be  noted  that  accurate  data  could  be  ob¬ 
tained  only  at  pressures  within  95  percent  of  collapse. 

These  results  are  of  some  interest  as  they  relate  to  the  problem  of 
"snap -through"  in  thin  cylindrical  shells.  The  Southwell  method  is  based 


iL 

In  practice  the  Southwell  method,  when  applicable,  predicts  a  buckling 
load  slightly  higher  than  that  actually  attained.  This  is  to  be  expected 
since  the  method  gives  the  elastic  buckling  load  for  a  geometrically  per¬ 
fect  structure. 
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Maximum  Stresses  Measured  during  Final  Run 
{Bay  4  at  110  deg,  630  psi) 


x/ L 

Stress 

Orientation 

Stress  in  psi 

Measured 

Pulos- 

Salerno 

Theory 

1,000 

Long.  Interior 

- 

72,900 

0.943* 

Long.  Interior 

52, 100 

0,500 

Circ.  Exterior 

74, 100 

Elastic 

Limit  =  74,400  psi 

*Longitudinal  gage  located  as  close  to  frame  as  possible. 

on  small  deflection  theory  and  predicts  the  buckling  load  for  the  geomet¬ 
rically  perfect  structure.  Hence,  it  should  not  give  accurate  results  for 
Sdi  tftirUdtures  Where-  tile  buckling  strength  can  only  be  explained  on 


the  basis  of  large  deflection  theory.  In  such  cases, "snap-through"  occurs 
at  some  load  less  than  the  load  predicted  by  small  deflection  theory.  The 
Southwell  results,  therefore,  indicate  that  if  "snap-through"  took  place,  it 
did  not  significantly  affect  the  buckling  strength. 

It  is  also  apparent  that  the  failure,  if  not  initially  elastic,  was  so 
nearly  so  that  the  collapse  pressure  wac  not  seriously  affected  by  inelas¬ 
tic  behavior.  Had  this  not  been  the  case,  the  Southwell  plots  would  have 
been  nonlinear  and  the  elastic  buckling  pressure,  though  possibly  obtain¬ 
able  from  strains  in  the  elastic  region,  would  not  have  agreed  well  with 
the  experimental  collapse  pressure. 


In  this  regard,  it  is  useful  to  examine  the  stresses  determined  from 


strains  measured  in  critical  areas.  The  maximum  stresses  were  observed 


in  Bay  4  at  110  deg  and  are  listed  in  Table  10  along  with  the  corresponding 
values  from  the  theory  of  Fulos  and  Salerno.33  All  stresses  are  for  a  pres¬ 
sure  of  630  psi.  The  differences  between  the  measurements  and  the  calcu¬ 
lations  can  be  attributed  to  additional  bending  stresses  which  accompanied 
the  lobar  deformations .  It  will  be  seen  that  none  of  the  measured  stresses 
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Figure  17  -  Distribution  of 
Nonlinear  Circumferential 
Strain  for  Bay  3  at  630  psi 
along  90-deg  Generator 


Figure  18  -  Theoretical  and 
Experimental  Distribution 
of  Linear  Strain  for  Bay  3 
(Run  Z,  Final  Test) 


reached  the  elastic  limit.  However,  it  should  be  pointed  out  that  the  in¬ 
terior  longitudinal  stress  axactly  at  the  juncture  of  frame  and  shell  could 
not  be  measured  because  the  longitudinal  gages  were  necessarily  displaced 
slightly  (G.057L)  from  the  juncture  to  allow  for  the  lexigth  of  a  gage  ele¬ 
ment.  The  maximum  strews  measured  at  this  neighboring  point  and  the 
theoretical  values  for  this  region  (Table  10)  indicate  that  the  elastic  limit, 
may  have  been  exceeded  by  the  interior  longitudinal  stress  in  Bay  4  at  110 
deg  .immediately  adjacent  to  Frame  5.  Because  of  the  'high  stress  gradi¬ 
ent  in  this  region,  a  precise  estimate  of  this  stress  was  not.  possible.  In 
any  event,  it:  appears  that  whatever  yielding  may  have  occurred  was  very 


[Si  Si 


slight  and  highly  localized. 

The  circumferential  gag ea  in  Bay  3,  located  along  the  90-deg  gener¬ 
ator,  provide  «ome  information  regarding  the  longitudinal  profile  of  a  lobe. 
Figure  17  thows  a  nondime  ns  ioonal  plot  versus  x/Lr£  for  a  pres- 

sure  of  630  psi.  A*  explained  previously,  t*  is  the  nonlinear  or  lobar 
component  of  the  total  strain,  is  the  value  of  e*  at  midbay,  and  x  is 
the  distance  measured  from  the  center  of  Frame  3,  Also  tfhowa  is  a  aine- 


wave  distribution  such  as  would  exist  in  the  case  of  simple  support.  The 
departure  of  the  points  from  this  sine  wave  indicates  to  some 

extent  the  influence  of  the  rotational  restraint  on  the  buckling  configura¬ 
tion,  Undoubtedly,  experimental  error  is  responsible  for  some  of  the 
deviations.  Strains  measured  on  Frame  4  in  this  region  showed  some 
evidence  of  asymmetric  deformations,  but  the  strain  variations  were  not 
large  enough  for  the  determination  of  a  well-defined  pattern. 

Straine  measured  in  Bay  3  also  provided  an  opportunity  for  further 
evaluation  of  the  axi symmetric  stress  analysis  of  Polos  and  Salerno**  in 


the  linear  range.  Figure  18  compares  the  theoretical  distribution  with 
•  the  distribution  o£  strain  sensitivities 


measured  across  Bay  3  at  90  deg.  Agreement  between  theory  and  exper¬ 


iment,  it  will  be  seen,  was  generally  good. 


COMPARISON  WITH  FREVIOUS  TESTS 

Table  11  summarizes  the  test  results  of  Cylinder  BR-4B  and  the 
earlier  results  for  Cylinders  BR-4  and  BR-4A,  The  cylinders  were  geo¬ 
metrically  similar  but  differed  in  their  yield  strengths.  The  table  in¬ 
cludes  pressures  given  by  Formula  [92A]  of  Von  Sanden  and  Gunther,3* 
at  which  the  exterior  circumferential  stress  at  midbay  reaches  the  yield 
value.  Since  this  pressure  for  BR-4  and  BR-4A  was  less  than  the  elastic 
buckling  pressure  given  by  Equation  £  33 j,  it  ia  clear  that  both  of  these 
failures  were  initially  inelastic.  Use  of  much  higher  strength  steel  in  the 
case  of  BR-4B  resulted  in  a  yield  pressure  well  above  the  elastic  buck¬ 
ling  pressure  and  in  a  mode  of  collapse  which  appears  to  have  been  ini¬ 
tially  elastic. 


Table  11 


Experimental  Results  Compared  with  Previous  Tests 


r 

Cylinder  Number 

— 

BR-4 

BR--4A 

BR-4B 

Fabrication 

Welded 

Machined 

Machined 

ffy»  psi* 

50,000 

50,  000 

82, 500 

®  i 

5*2 

jj  J 

2  •  » 
•3  g  o. 

U(S 

Von  Sanden 
and  CKinther 

Formula  [92A],Z 

494 

494 

815 

Equation  [33] 

633(11)**  | 

Experiment 

390(10) 

550(10-11) 

633(11) 

*Based  on  an  offset  strain  of  0.002  . 

**Number  of  lobes  in  parantheses. 

* 

Dimensions  and  Young's  Modulus  (28.9  x  10°  psi) 
r :  ‘  as  sumed  identical  for  the  •  three  cylinder  rr  -  * 


SUMMARY  AND  CONCLUSIONS 

1.  The  test  Cylinder  BR-4B  collapsed  at  633  psi  in  an  asymmetric 
shell  buckling  mode  characterized  by  11  circumferential  lobes.  On  the 
basis  of  strain;  measurements,  it  is  concluded  that  the  collapse  was  - 
initiated  by  elastic  buckling. 

2.  The  Southwell  method  can  be  an  effective  means  of  determining 
nondeatructivcly  the  elastic  buckling  strength  of  short  shells.  It  is  likely 
to  be  accurate  only  where  imperfections  and  residual  stresses  are  small, 
as  in  the  case  of  machined  cylinders,  and  only  where  it  is  possible  to  ap¬ 
proach  very  closely  the  elastic  buckling  pressure.  Hence,  its  practical¬ 
ity  aa  a  nondestructive  technique  may  be  rather  limited  for  this  mode  of 
buckling . 
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PART  HI  -  EVALUATION  OF  BUCKLING  THEORY 


OH 
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COMPARISON  OF  THEORY  WITH  EXPERIMENT 


According  to  the  theory  developed  in  this  report  (Equation  [33]), 
buckling  for  Cylinder  BR-4B  occurs  at  633  psi  in  a  configuration  of  II 
circumferential  lobes,  which  is  in  exact  agreement  with  the  test;  results. 
Such  agreement  is,  of  course,  better  than  one  has  any  right  to  expect, 
even  for  the  most  rigorous  theory  imaginable.  The  present  theory  makes 
use  of  several  approximations;  moreover,  its  accuracy  depends  on  the  uni¬ 
formity  of  the  test  cylinder  and  the  accuracy  with  which  its  properties  can 
be  determined.  The  shell  thickness,  for  example,  must  necessarily  be 
represented  by  an  average  of  many  measurements,  all  subject  to  error. 
The  determination  of  Young's  modulus  required  a  separate  test  which  uti¬ 
lized  the  response  of  a  ring  to  diametral  loading.  While  this  method  is 
believed  reliable,  it  is  still  subject  to  experimental  error,'  and  no  claim 


is  made  that  it  is  the  best  that  could  have  been  used.  Other  techniques, 
for  example,  theuse  bf Ispllsajr 'jtffih' ^|es“ W  rneaoure  tbe  response  of^’^*-'  • 
a  specimen  under  direct  stress,  are  highly  regarded.  Because  perfect 
isotropy  is  never  achieved,  methods  which  employ  different  conditions  of 
stress  cannot  be  expected  to  yield  identical  results. 

In  view  of  these  uncertainties,  the  degree  of  correlation  of  theory 
with  experiment  can  be  regarded  with  some  suspicion,  but  the  results  are, 
nevertheless,  substantial  evidence  of  the  validity  of  the  theory. 

The  test  of  BR-4B  is  considered  successful  as  a  demonstration  of 


the  phenomenon  of  elastic  buckling,  but.  it  is  not  an  ideal  example  of  the 
influence  of  stiffening  ringn  on  shell  buckling  strength.  The  buckling  pres¬ 
sure  was  only  slightly  greater  than  that  (600  poi*)  given  by  the  theory  of 
Von  Mises  (Equation  [46])  fur  the  case  of  simple  supports.  This  is  a 


*This  figure  is  obtained  on  the  basis  that  buckling  is  confined  to  the 
unsupported  length  (L)  of  shell  plating.  If  the  full  frame  spacing  (  Lr) 
is  used,  the  resulting  pressure  is  560  psi,  as  indicated  in  Table  4. 
Still  less  is  the  pressure  (550  psi)  given  by  Equation  [42],  which  in¬ 
cludes  the  effect  of  the  boundary  conditions  on  the  initial  deflections. 
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natural  consequence  of  the  large  value  of  (3  (6.  80)  used  in  an  effort  to 
obtain  an  elastic  failure.  Nevertheless,  the  test  did  demonstrate  that 


the  Von  Mises  pressure  can  be  exceeded. 

Fortunately,  since  the  test  of  BR-4B,  other  results  have  become 
available  that  provide  better  data  regarding  the  influence  of  stiffening 
rings.  These  have  been  selected  from  tests  of  a  number  of  machined 
steel  cylinders,  8  in.  in  diameter,  most  of  which  were  recently  reported. 
The  cylinders  in  many  cases  had  more  typical  bays  than  did  BR-4B,  but 
there  is  no  indication  that  the  extra  length  had  any  influence  on  collapse 
other  than  to  isolate  the  central  bays  from  end  effects.  In  all  cases, 
failures  were  of  the  interframe  variety,  but  in  only  a  few  was  there  evi¬ 
dence  of  elastic  buckling.  Since  none  of  them  was  instrumented  for  the 
specific  purpose  of  datecting  an  elastic  failure,  the  cases  cited  here  were 
selected  by  comparing  pc,  the  critical  buckling  pressure  according  to 
Equation  [33],  with  an  estimated  yield  pressure  py,  given  by  Von  Sanden 
and  Giinther  Formula  [92A].JZ  The  pressure  py  is  that  at  which  the  ex¬ 
terior  circumferential  stress  at  midbay  reaches  the  yield  value  defined 

at  an  offset  strain  of  0,002.  When  p  exceeds  p  ,  the  possibility  of  elas- 

Y  ,  c 

tic  buckling  exists.  If  p  exceeds  p  ,  the  failure  must  be  inelastic.  On 

c  y 

this  basis,  three  of  the  cylinders  were  in  the  elastic  category.  .Their  prop¬ 
erties  and  those  of  BR-4B  are  listed  in  Table  12.  Since  modulus  meas¬ 
urements  for  the  small  cylinders  were  not  available,  a  nominal  value  of 
30  x  10^  psi  was  assumed.  The  collapse  pressures,  both  experimental 
and  theoretical,  for  the  four  cylinders  appear  in  Table  13, 

The  influence  of  the  frames  is  clearly  seen  by  comparing  the  Von 
Mises  pressures  in  the  table  with  the  experimental  results.  The  additional 
strength  provided  by  the  frames  (as  much  as  22  percent)  results  prima¬ 
rily  from  the  fact  that,  for  each  of  the  small  cylinders,  (3  was  considerably 
less  than  for  BR-4B.  Evidently  the  present  theory  adequately  accounts 
for  this  effect.  In  each  case,  the  pressure  given  by  Equation  [33]  agrees 
well  with  the  experimental  pressure.  These  results  are  also  presented 
graphically  in  Figure  19.  The  abscissa  is  pc/p  and  the  ordinate  is  the 
ratio  oi  p,  ,  the  experimental  failure  pressure,  to  p^..  The  4f»-deg  line 
drawn  on  the  figure  represents  perfect  agreement  between  theory  and 


experiment. 


U-  T>C/Py 


is  less  than  unity  (indicated  by  the  dotted  line), 
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Table  12 


elastic  buckling  ia  expected.  It  will 
be  seen  that  the  four  cases  of  Table 
13  all  fall  within  this  region  and  the 
plotted  points  all  lie  close  to  the  45- 
deg  line.  Some  of  the  other  tests 
reported  in  Reference  39  are  also 
plotted  in  the  figure  along  with  the 
test  of  BR-4A.  In  all  of  these  cases, 
p./p.r  exceeds  unity,  so  that  the 
points  are  outside  the  elastic  region 
and,  as  expected,  all  lie  below  the 
45-dog  line. 

These  results  indicate  that  the 
performances  of  the  four  cylinders 
-have  been  satisfactorily  explained 
using  small-deflection  theory,  and  are  consistent  with  earlier  remarks 
regarding  the  application  of  the  Southwell  method  for  the  case  of  BR-41V 
w  This  doe  s  not  nec  e  s  »ar  ily  Indicate  that  the  Observed  failure s  we r  e  heti  lii  $ 
fact,  precipitated  by  a  "snap- through"  mode  of  instability.  That  this  phe¬ 
nomenon  can  occur  in  cylindrical  shells  under  hydrostatic  pressure  has 
been  demonstrated  experimentally  by  Kirstein  and  Wenk,21  and  theoretical¬ 
ly  by  the  large  deflection  analyses  of  Donnell,17*18  Langhaar  and  Boresi,17 
and  Kempner  and  Crouzet-Pascal,20  among  others.  However,  these  studies 
have  also  shown  that  the  pressure  at  which  "snap-through"  takes  place  is 
not  much  Iobb  than  the  pressure  given  by  Bmall-deflection  theory  even  in 
the  presence  of  imperfections.*  For  machined  shells,  this  pressure  dif¬ 
ference  might  be  too  small  to  be  detectable. 

Table  13  also  shows  that  according  to  the  solution  of  Von  Sanden  and 
Tolke,12  the  pressures  reflect  some  strengthening  due  to  the  presence  of 
the  frames,  but  still  fall  short  of  experiment.  This  is  to  be  expected  since 


Figure  19  -  Experimental  and 
Theoretical  Results 


*This  is  in  marked  contrast  to  the  behavior  of  cylindrical  shells  under 
axial  compression  or  torsion  where  the  pressure  difference  can  be  very 
large. 
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the  analysis  considers  the  effect  of  the  frames  on  the  prebuckling  defor¬ 
mations  but  fails  to  account  for  the  restraint  which  they  provide  during 
buckling.  On  the  other  hand*  the  Nash  solution  for  a  obeli  with  lull  fixity^® 
trreatlv  overestimates  the  strength  actually  realise,'}  in  each  rase,  as  the 
table  indicates.  This  should  be  interpreted  not  as  a  refutation  of  the  solu¬ 
tion,  but  rather  as  evidence  of  its  inapplicability  for  the  cases  considered. 
It  is  possible,  as  has  already  been  indicated,  that  the  analysis  may  be 
quite  accurate  in  cases  where  complete  fixity  is  actually  attained. 

Finally,  it  should  be  said  that  the  experimental  results  reported 
here,  though  meager,  are  all  of  which  the  author  is  aware.  Unfortunately, 
there  are  no  data  for  internally  framed  cylinders,  and  none  are  apt  to  be 
forthcoming  because  of  the  machining  problems  presented  by  internal 
frames. 

While  no  future  studies  in  the  area  of  elastic  shell  instability'  are 
contemplated,  should  additional  experimental  data  become  available,  they 
will  be  examined  with  interest. 

■  >'•  . . .  ■  — .  -v  v;. w7;., 

CONCLUSIONS 

1.  The  small  deflection  theory  of  this  report  predicts  with  accuracy 
the  elastic  buckling  of  cylindrical  shells  supported  by  closely  spaced  ring 
stiffeners,  at  least  where  stiffeners  are  external  and  p  is  greater  than 
4.0.  Because  of  the  restraint  provided  by  the  rings,  the  buckling  strength 
can  be  considerably  greater  than  that  expected  for  a  shell  with  simple 
support. 

2.  The  solution  of  Von  Sanden  and  Tdlke  is  conservative  in  all  cases 
because  it  accounts  only  partially  for  the  strengthening  influence  of  the 
stiffeners. 

3.  The  fact  th&t  the  Nash  analysis  for  a  shell  with  clamped  edges 
greatly  overestimates  the  buckling  strength  of  ring- supported  shells  does 
not  necessarily  indicate  the  use  of  an  unrealistic  buckling  shape  in  the 
analysis.  It  is  more  likely  that  the  disparity  exists  chiefly  because  frames 
of  practical  size  do  not  provide  restraint  comparable  to  complete  fixity. 

4.  Since  the  Von  Mises  analysis  is  based  on  the  assumptions  that  the 
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shell  is  unsupported  prior  to  buckling  and  simply  supported  during  buck¬ 
ling,  it  can  be  highly  inaccurate  in  predicting  buckling  pressures  for  i ing- 
stiffened  shells.  Nevertheless,  it  is  probably  the  most  practical  means 
for  estimating  elastic  buckling  strength  because  it  is  always  conservative 
and  can  be  represented  in  a  very  simple  form. 
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APPENDIX  A 


STRAIN  ENERGY  OF  THE  SHELL 


In  the  development  which  follows,  basic  strain-displacement  rela¬ 
tions  are  obtained  with  accuracy  preserved  through  quadratic  terms  in 
the  displacements.  The  strain  energy  integral  of  the  shell  is  then  formu¬ 
lated  with  displacement  terms  maintained  through  the  third  order. 

We  consider  an  element  situated  within  the  thickness  of  the  shell 
at  a  radius  r  and  having  axial  and  angular  coordinates  a  and  6.  The 
element  undergoes  displacements  u' ,  v\  and  w1  in  the  axial,  tangential, 
and  radial  directions,  respectively,  with  w'  taken  as  positive  Outward. 

We  will  ss suit  j  that  the  strain  in  the  radial  direction  is  negligible  so 
that  only  deformations  within  the  plane  of  the  element  need  be  considered. 
Novozhilov40  describes  these  deformations  in  terms  of  strain  components 
which  will  be  designated  ,  and  «X0  for  cylindrical  coordinates. 

They  are  related  to  the  displacements  as  follows: 


XX 


=  u 


x+It(«X2+  +  *0 


,.2 


,.2* 

'x 


vb 


+  w' 


“90  - 


UA  /Ufl  V  /VA  +  W'\  /W0-V'\ 

‘*8  =  vx+T  +  ui(T)  +  vi(~)+ W*(^) 


where  the  subscripts  on  u",  v',  and  w'  indicate  differentiation.  In  order 
to  calculate  strain  energy,  it  is  necessary  to  consider  strains  of  line  ele¬ 
ments  that  lie  along  mutually  perpendicular  axes  -  in  this  case  the  x  and 
0  coordinate  lines.  These,  strains,  according  to  Novozhilov,40  can  be  ex¬ 
pressed  in  terms  of  the  strain  components  by 


ex  "  ^1  +  2  x  -  1 
eQ  =  Vl  +  2  (an  ~  1 


[67] 


sin"'  t 


[67] 

continued 


ex6 


x0 _ 

^0“+  2*^(1  +  2<0e) 


where  ex  and  eg  arc  unit  strains  in  the  x  and  8  directions  and  ex0  is 
the  shear  strain  expressing  the  change  in  the  angle  between  the  coordinate 
lines.  By  combining  Equations  [66]  and  [67]  and  discarding  all  displace- 
ment  terms  of  order  higher  than  quadratic,  one  obtains 


e 


x 


“it  +  J  <vic2  + 


v^) 


e8 


ve  +  w 


[68] 


®X0 


ue  ,  ,/wa-v\  ,  .  ue/ve+w'\ 

_  +  wi  f—_)  -  »ivi  -  —  (—7—) 


uk  ,vk  +  w\ 


It  may  be  surprising  to  note  that  retention  of  only  linear  terms  in  the 
strain  components,  as  was  apparently  done  by  Langhaar  and  Boreai,19 
leads  to  a  slightly  more  complicated  result.* 

To  relate  the  displacements  u’,  v',  w'  at  any  point  to  the  displace¬ 
ments  u,  v,  w  at  the  middle  surface  of  the  shell,  it  is  useful  to  employ 
relatione  developed  by  Langhaar  and  Boreai. 


u*  =  u  I 

V1  V  +  7.U  [69] 


w'  =  w  +  zY 


In  that  case,  the  ati-aiuo  are  identical  with  the  corresponding  strain 
components.  However,  the  difference  between  Equations  [66]  and  [68] 
ia  probably  of  little  practical  significance. 
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where 


'w0 "v\ 

£  -  -  w  !  v„  I - 1  i  u  v. 

s  X  R  J  XX 

(wp  -  V .  UA  w  >vfi+wwwA-v. 

-V)  + V+(V)(-V) 


V  = 


*4  /wo'v\2 


X  / 

2  "V 


R  / 


where  R  is  the  radius  to  the  middle  surface  of  the  shell  and  z  is  the 
thickness  coordinate  measured  positive  outward  from  the  middle  surface 
of  the  shell.  Equations  [69]  are  based  on  the  assumption  of  Kireboff  that 
normals  to  the  undeformed  middle  surface  remain  normal,  straight,  and 

unextended  after  deformation.  By  combining  Equations  [6ff]  and  [69] 
after  suitable  differentiation,  the  resulting  equations  at#- 


ex  *  Ujj  '  * +  i  {v2  +  w2)  +  *(Nvxx  +  uxxwx  +  + 


z2N2 


•  *"e  i  /«!  rZ\  z  /"ee"* 

-  — - —  +  ~  (-=■  4  N^1  +  —  ( -  +  MftN  + 

R  2  \n2  /  R  \  R  6 


^  0  y  .1 

ex6  =  vx  +  ~  -R  (RNx+  wx6)  +  wxN  -uxvx-— '  +iH2ux9wx 


R  R 


+  uGwxx+  uxw*e  4  RMxN+  RMNx  +  N8vx+Nvx6  +  wxN 
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-*■  RuvN  + - —  +  w  .  M  ~  (RZ  w„„N„  +  w  e  Ng) 
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where 
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v0  +  w 
R 


N 


WA  "v 
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In  obtaining  Equations  [70],  it  was  assumed  that  r  could  be  replaced  by 
It  with  negligible  loss  of  accuracy.  A  further  simplification  can  be  made 
in  exQ  by  neglecting  z/R  in  comparison  to  unity .  Then 

cx0  *  vx+  r-“r  <RNx+wx8)  +  wxN*'uxvx“—  [71) 


+  7  (2ux0wx  +  ^  wxx  +  uxwx0  +  RMxN  +  RM  Nx  +  Nevx 

+  N  vx0  +  R Nx  +  ^  *  wx9  ^ R  wxx Nx  +  wx6 

In  formulating  the  expression  for  the  strain  energy  of  the  shell,  it 
will  be  assumed  that  the  stress  in  the  radial  direction  is  zero.  According 
to  Langhaar  and  Boresi,19  the  shell  energy  Ua  can  then  be  expressed  by 


,2v 


U8  = 
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tfy  S  I  I h  ['*  +  «e  +  2  v  e*  '0  +  ('T1)  ] R  dxde  dB 

X'  A  V*\ 


’ n  -  — 
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where  E  is  Young's  modulus  and  h  is  the  shell  thickness.  After  combin¬ 
ing  Equations  [70]  and  [71]  with  Equation  [72],  performing  the  integra¬ 
tion  on  %  and  discarding  terms  higher  than  the  third  order,  one  obtains: 
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u2 
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APPENDIX  B 


STRAIN  ENERGY  OF  RINGS 


The  strain  energy  of  a  ring  will  be  considered  as  the  sum  of  three 
quantities:  Ujp  ,  the  energy  due  to  deformations  in  the  plane  of  curva¬ 
ture;  Uop  ,  the  er  'rgy  of  deformations  perpendicular  to  the  plane  of 
curvature;  and  U\j,qR  ,  the  energy  due  to  torsion  about  an  axis  parallel 
to  the  0 -coordinate  line.  Thus,  the  strain  energy  of  a  ring  is 


uf  =  uIp  +  uop  +  uTOR 


[74] 


In  this  equation,  the  energy  associated  with  restricted  warping  as  well  as 
that  due  to  extensions  of  the  ring  in  the  a  and  x  directions  has  been  neg¬ 
lected. 


Ujp  is  simply  the  energy  associated  with  the  circumferential  strain, 


or 
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[75] 


where  Af  is  the  cross-sectional  area  of  a  ring  situated  at  a  distance  iLj 
from  the  origin.  From  Appendix  A,  Equation  [68]  is 


e0 


ve  +  w'  if/  "e  \2 
R  +  e  2  [vR  +  e/  \  R+e  /  J 


[76] 


when  the  approximation  is  made  that  r  can  be  replaced  by  R+e.  Refer¬ 
ring  to  Figure  3,  e  is  the  distance  between  the  middle  surface  of  the  shell 
and  the  center  of  gravity  of  the  ring  cross  section.  When  substitutions 
from  Equation  [69]  of  Appendix  A  are  made,  Equation  [76]  becomes 


After  this  expression  is  substituted  in  Equation  [75]  and  terms  higher 
than  the  third  order  are  discarded,  one  obtains 
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Before  integrating  over  the  area  A^ ,  it  is-  helpful  to  make  use  of  the  fol- 
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is  the  moment  of  inertia  about  an  axis  normal  to  the  area  A £  at  its 
center  of  gravity.  Ujp  can  then  be  written: 
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According  to  fundamental  beam  theory,  the  out-of-plane  bending  energy 
is 

El. 
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[81] 


where  Izq  is  the  moment  of  inertia  of  the  ring  about  an  axis  parallel  to 
the  z-coordinate  through  its  center  of  gravity  and,  according  to  Reference 
34, 


X  = 
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182] 


Uq  is  the  axial  displacement  of  the  centroid  of  the  ring;  £  is  the  rotation 
of  the  ring  in  the  x-z  plane  as  defined  in  Equation  [69]  of  Appendix  A. 
Because  each  ring  is  assumed  to  undergo  no  strain  in  the  axial  direction, 
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With  x  replaced  by  e  in  Equation  [69]  of  Appendix  A,  it  follows  that 
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Combining  Equations  [81]  and  [84]  and  discarding  terms  higher  than  the 

third  order,  one  obtains  • 
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According  to  Reference  34,  the  torsional  energy  is  expressed  by 


U  =  ek 
TOR  4  (  1+v) 


$>  (R  +  e)  dG 


where  $>  is  the  angle  of  twist  per  unit  of  circumference  and  K  is  the  tor¬ 
sion  constant.  For  example,  in  the  case  of  a  rectangular  ring  where  the 

i  3 

depth  d  is  much  greater  than  the  width  b,  K  is  given  by  f  db  .  When  an 
open  section  is  comprised  of  several  narrow  rectangles,  K  for  the  section 
is  given  approximately  by  tire  sum  of  the  values  K  for  each  part.4*  Refer¬ 
ence  42  gives  values  of  K  for  rectangles  where  d  does  not  greatly  exceed 
b.  Reference  34  defines  the  angle  of  twist  as 
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After  appropriate  substitution  this  becomes 


(R+e) 


-  [R(-wrt  +  vxeN  +  vxNe)  -  “e  J  x  = 


and  the  torsional  energy  is 
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with  terms  above  third  order  discarded. 
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APPENDIX  C 


POTENTIAL  OF  THE  EXTERNAL  LOADS 


The  wui k  June  by  Lhe  external  ioade  is  defined  as  tfte  product  ot  tbe 
pressure  and  the  change  in  volun.e  undergone  by  the  cylinder  during  defor¬ 
mation,  or 


"W  =  p(V'-  V)  190] 

where  V  and  V'  are  the  volumes  before  and  after  deformation,  respec¬ 
tively.  The  volume  bounded  by  an  element  of  the  deformed  shell  is  defined 

dV'  =  dA'df  [91] 

where  A'  is  the  area  enclosed  by  the  median  surface  of  the  shell  in  the 
R  -  G  plane,  and  df '  is  the  length  of  the  element  measured  in  the  ^direction. 


Figure  20  -  Deflected  Element  of  Shell 
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The  coordinates  can  be  expressed  as  follows: 
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df  '  =  (1  +  u-Jdx  . 

When  quantities  above  the  second  order  are  neglected, 

dA*  =;  [(R  +  w)2  +  vfl(R  +  w)  +  v(v+  w-)] 


[94] 


and 
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But,  since  J  [  V0  *  ~r  (v$  w  *"  vwq)J  *3®  vanishes  because  of  c 
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APPENDIX  D 


if 

St 


INFINITE  SERIES  FORMULATION  OF  THE  INITIAL  DEFLECTIONS 


»<r- 


i 

* 

I 


i  he  axi»ymiiiei,jriC  daouc  defer iiiatiOnS  of  a  ring  -  stiffened  cylindri 
cal  shell  under  external  pressure  were  first  analyzed  by  Von  Sanden  and 
Gunther, 32  This  work  was  later  extended  by  Pulos  and  Salerno33  to  include 
the  nonlinear  or  "beam-column"  effect  of  the  axial  load.  To  express  the 
deflections  in  the  form  of  a  trigonometric  scries  it  would,  therefore,  be 
possible  to  perform  a  straightforward  Fourier  expansion  of  either  of  these 
closed  form  solutions.  An  alternative  approach,  however,  is  favored  here' 
since.it  leads  quite  simply  to  the  desired  result  directly  from  the  equa¬ 
tions  of  equilibrium  with  no  dependence  upon  prior  solutions. 

The  ajd symmetric  deformations  of  the  cylinder  are  described  as 
follows; 

'• _  _ _  _  U  «  u(x) 


W  -  W  { x) 

with  the  added  condition  that  wx  vanishes  at  a  frame.  Referring  to  Equa¬ 
tion  [2]  then,  we  have  for  the  initial  total  potential. 
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Products  of  hJ  with  third-order  terms  in  the  displacements  have  been  neg¬ 
lected  in  this  equation.  Because  deformations  in  each  bay  are  identical, 
it  is  only  necessary  to  carry  out  the  integration  on  x  from  0  to  Lf  . 
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...  ,^4— 


Integration  with  respect  to  8  can,  of  course,  be  eliminated  because  of 
axial  symmetry. 


If  we  consider  the  dependent  variables  to  be  uv  and  —  ,  then  one  of 

x  R  _ 

the  equations  which  results  from  equating  the  first  variation  o£  U-j.  to 
zero  is 
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199] 


By  making  this  substitution  in  Equation  [98],  retaining  only  linear  and 
quadratic  terms  in  w  and  its  derivatives,  and  neglecting  products  of  p‘ 
with  w,  one  obtains 
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w  is  now  taken  to  be 
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where  X_  is 

Lj 

and  the  integration  is  performed 


When  this  series  is  substituted  into  Equation  [100] 
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The  equations  of  equilibrium  are  obtained  by  minimising  with  respect  to 
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where  ft*  is  ——— ~ — - —  .  Subtracting  Equation  [104]  from  Equation  [103] 
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when  — —  is  negligible  compared  with  unity  and  2X^  »  1  -  2v.  Summing 
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the  Wj'>j  we  obtain: 
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w  ill  then  determined  by  combining  Equations  (103]  and  [106] 
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am  can  Put  *n  a  TOOr®  convenient  form  as  follows:  If  an  unstiffened  cylin- 


/2Xmx\ 

d«r  were  to  buckle  in  the  axi  symmetric  configuration  w  a  am  cos  I  — — — — 1 , 
the  corresponding  buckling  pressure  pm  would  be  given  by 


'  •--•"•• '  . . “  “■  -  ~p  "b  i  riV., 

*rn  #— 
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This  equation  a! no  applies  in  the  case  of  a  long  cylinder  under  end  thrust 
where  the  length  of  a  complete  axial  buckling  wave  is  equal  to  L^/m,  If 
Equations  [107]  and  [108]  are  combined,  we  find  that 
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From  Equations  [105]  and  [107],  we  have 


w. 


w 


PR  /-> 

-  a  ■  — —  ( 2  -  v  /  r>, 
J  2Eh 


pR/'  . 
- - (Z  -  v) 


[■-’ 


Zr) 


am  COfl 


m~  1 


/  ^  x  \  1 


f  1 10] 


92 


w  = 


(2-v) 


n 


[no] 

continued 


since  aQ  equals  unity.  The  combination  o£  Equations  [110]  and  [99] 
gives 


1 1 

“x 


r, 

2Eh  L* 


,  /  •> 
\  ** 


m 


f  2*mx 

\  R 


m=  -«o 


if  quadratic  terms  in  p  are  neglected. 

One  consequence  of  Equation  [110]  is  a  condition  for  axisymmetric 
buckling.  It  can  be  seen  that  w  will  increase  without  limit  should  1/q 
approach  aero.  The  governing  condition  then  is 


So  long  as  A^  >  0,  this  condition  can  be  met  only  when  p  >  pm  for  one  or 
more  values  of  m.  If  Aj  =  0,  Equation  (112]  is  satisfied  when  p  »  pm> 
This  is  to  be  expected  since  Equation  [108]  holds  only  for  buckling  of  an 
unstiffened  cylinder.*  In  that  case  the  expression  in  Equation  [110'j  for 
w  becomes  indeterminate. 

For  purposes  of  comparison,  the  closed-form  solution  of  Pul  os  and 
Salerno33  for  w  can  be  written  as  follows: 


pR 
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*It  should  be  emphasized- that  the  buckling  configurations  mentioned 
here  are  all  of  a  single  class  characterized  bv  symmetry  about  a  frame 
and  need  not  correspond  to  the  minimum  buckling  pressure.  Configura¬ 
tions  of  the  nonsymmetric  class .  permitting  frame  rotation,  can  result 
in  lower  pressures,  depending  on  the  geometry  of  the  cylinder. 
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where 


f(x)  =  [jl^  [  cosh  x  sinK^  (L  -  x)  +  cosh  Kj  (L  -  x)  sin  K2  x  ] 
+  n2  [  sinh  Kj  x  cos  K-2  (i-  -  x)  +  sinh  Kj  (L  -  x)  cos  K2  x  ] 
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where  L  ia  the  inner,  or  unsupported,  frame  spacing  and  b  is  the  faying 
width  of  a  frame.  If  b  ia  neglected,  L  equals  Lf  and  liquations  [110] 
and  [113]  should  be  equivalent.  Equating  these  two  expressions  at  x  =  0, 
we  find  that 
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Combining  Equations  [113]  and  [115]  yields  tjjie  additional  result: 
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The  principal  reason  for  introducing  the  Pulos  and  Salerno  solution 
was  to  provide  the  alternative  means  (Equation  [  114])  wliemneeded  for 

CO 

computing  am.  Graphs  of  the  function  G  can  be  found  in  Reference 

me  “00 

33  where  it  is  designated  Fj  for  p  <  p*  and  Fg  for  p  >  p*.  In  cases 
where  p«  p*,  it  follows  from  Equation  [108]  that  p«  pm«  Then  Equa¬ 
tion  [114]  reduces  to 
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1  ,  -  ^ 

For  purposes  of  buckling  calculations*  tills  latter  equation  is  usually  suf¬ 
ficient.  Furthermore*  it  is  often  as  simple  to  work  with  t)*e  series  form 
since*  except  for  large  values  of  p,  convergence  is  quite  rapid.t  Usually 
only  three  or  four  terms  will  provide  good  accuracy. 


tu  sing  a  compttrison  tost  with  the  convergent  series  23 
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one 


m- 1  m' 


can  easily  show  that  the  series  (Equation  [  117])  converges  for  all  finite 
values  of  (3. 
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APPENDIX  E 


STRESS  FUNCTION 


The  energy  expression  in  Equation  [  5J  is  extremely  unwieldy  as  it 
stands.  It  can  be  greatly  simplified  if  u  and  v  can  be  somehow  eliminated 
To  do  this,  use  will  be-  made  of  a  stress  function  following  the  method  of 
Donnell,43 

The  strain-displacement  relations  with  only  linear  terms  retained 
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But,  from  the  fundamental  definitions  of  plane  strain 
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where  ax,  <Jq  and  <7xg  are,  respectively,  the  axial,  circumferential, 
and  shear  stresses.  Combining  Equations  [119]  and  [120]  we  find 
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The  stress  function  F  is  defined  as  follows: 
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Equation  [121]  now  becomes 
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Once  w  has  been  specified,  F  is  determined  from  Equation  [123],  and 
u  and  v  are  obtained  from 
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APPENDIX  F 


DETERMINATION  OF  THE  MODULUS  OF  ELASTICITY 
FOR  CYLINDER  BR-4E 


Young's  modulus  for  the  test  cylinder  was  determined  by  measuring 
the  deflections  under  concentrated  loading  of  a  ring  taken  from  one  end  of 
the  cylinder.  Although  not  in  common  practice,  this  method  has  certain 
advantages.  One  is  that  relatively  large  elastic  deflections  can  be  pro¬ 
duced  which  are  easier  to  measure  accurately  than,  for  example,  strains 
in  an  axially  stressed  rod.  Another  advantage  is  that  a  larger,  and  there¬ 
fore  more  representative,  specimen  of  the  cylinder  can  be  used. 

After  the  test  of  BR-4B  had  been  completed,  the  bulkhead  ring  near¬ 
est  the  damaged  portion  (Figure  13)  was  cut  from  the  cylinder  and  ma¬ 
chined  to  the  following  dimensions: 


Mean  radius  (R)  --  8. 196  in. 


The  ring  was  placed  vertically  between  the  heads  of  a  universal  testing 
machine  and  was  set  in  position  by  an  initial  compressive  load  of  50  lb. 
The  load  was  then  increased  in  50-lb  increments  up  to  350  lb  yhile  the 
shortening  of  the  vertical  diameter  of  the  ring  was  measured  with  a  dial 
indicator.  The  ring  was  unloaded  in  the  same  manner.  Two  such  runs 
were  made  and,  in  each  case,  measurements  during  loading  and'  unload¬ 
ing  agreed  within  0.  002  in.  The  results  of  the  test  are  shown  graphically 


in  Fitmre  7. 1 

Hit 


According  to  Reference  44,  for  a  ring  of  rectangular  cross  section, 
the  change  A  in  the  diameter  produced  by  the  load  P  can  be  expressed 
as  follows; 


A  = 


PR^ 

Ebdc 


[125] 
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Tzryi ir-*- 


Otene*  m  Diom«t«r(A)  in  !nel*t  *  C'1 


Figure  21  -  Determination  of  Young^s  Modulus  for  BR-4B 
(Response  of  Ring  to  Diametral  Loading) 

The  equation  is  applicable  to  thick  rings  and  accounts  for  direct  and 
shear  stresses  as  well  as  bending  effects.  For  the  ring  in  question,  the 
result  is 

E  =  2. 88  x  104  ~  [  126] 

—  — - 

From  the  slopes  of  the  plots  of  Figure  21,  the  values  for  E  in  pounds  per 
square  inch  were  28.8  x  10^  and  29.0  x  10^.  An  average  value  of  28.9  xl0~ 
psi  was  thereby  determined. 
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